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Abstract 

We prove upper and lower bounds on the free energy of the Sherrington-Kirkpatrick model with multidi- 
mensional spins in terms of variational inequalities. The bounds are based on a multidimensional extension 
of the Parisi functional. We generalise and unify the comparison scheme of Aizenman, Sims and Starr and 
the one of Guerra involving the GREM-inspired processes and Ruelle's probability cascades. For this 
purpose, an abstract quenched large deviations principle of the Gartner-Ellis type is obtained. We derive 
Talagrand's representation of Guerra's remainder term for the Sherrington-Kirkpatrick model with multi- 
dimensional spins. The derivation is based on well-known properties of Ruelle's probability cascades and 
the Bolthausen-Sznitman coalescent. We study the properties of the multidimensional Parisi functional by 
establishing a link with a certain class of semi-linear partial differential equations. We embed the problem 
of strict convexity of the Parisi functional in a more general setting and prove the convexity in some par- 
ticular cases which shed some light on the original convexity problem of Talagrand. Finally, we prove the 
Parisi formula for the local free energy in the case of multidimensional Gaussian a priori distribution of 
spins using Talagrand's methodology of a priori estimates. 
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1. Introduction 

The Sherrington-Kirkpatrick (SK) model of a mean-field spin-glass has long been one of the most 
enigmatic models of statistical mechanics. The recent rigorous proof of the celebrated Parisi formula for 
its free energy, due to Talagrand [f30||, based on the ingenious interpolation schemes of Guerra jl7[ | and 
Aizenman, Sims, and Starr [[[J constitutes one of the major recent achievements of probability theory. 
Recently, these results have been generalised to spherical SK-models [^9[ and to models with spins taking 
values in a bounded subset of R rf2Tj]. 
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In this paper, we are mainly concerned with the question of the validity of the Parisi formula in the case 
where spins take values in a cZ-dimensional Riemannian manifold. We address the issue of extending the 
approach of Aizenman, Sims and Starr, and the one of Guerra to the multidimensional spins. We study the 
properties of the multidimensional Parisi functional. Motivated by a problem posed by plj], we show the 
strict convexity of the local Parisi functional in some cases. 

We partially extend Talagrand's methodology of estimating the remainder term to the multidimensional 
setting. In the case of the multidimensional Gaussian a priori distribution of spins we prove the validity of 
the Parisi formula in the low temperature regime. 

Definition of the model. Let E C R d and denote Eat = T. N . We define a family of Gaussian processes 
X = {X(a)}aeL N as follows 

X(a) =X N (a) = i £ gij(oi,Oj), d-D 
i'J=i 

where the interaction matrix G = {gij}f,- = i consists of i.i.d. standard normal random variables and, for 
x,y G R d , (x,y) = Y,i=i x uyu is the standard Euclidean scalar product. In what follows all random variables 
and processes are assumed to be centred. We shall call Hn{o) = —\/NXn(g) a random Hamiltonian of 
our model. 

Throughout the paper, we assume that we are given a large enough probability space (n,J?,P) such 
that all random variables under consideration are defined on it. Without further notice, we shall assume 
that all Gaussian random variables (vectors and processes) are centred. 

We shall be interested mainly in the free energy 

p N (P) = Uog^x V ^VNX(a))dii m (a), (1.2) 
where j3 > is the inverse temperature and ji € ./#f(E) is some arbitrary (not necessarily uniform or 



discrete) finite a priori measure. We assume that the a priori measure fi is such that (1.2) is finite. We shall 



be interested in proving bounds on the thermodynamic limits of these quantities, e.g., on 

p(P)=Kmj N (P). (1.3) 

Remark 1.1. Note that there is no need to include the additional external field terms into the Hamiltonian 
(JTTTJ), since they could be absorbed into the a priori measure ji. 

Mean-field spin-glass models (see, e.g., [0]) with multidimensional (Heisenberg) spins were considered 
in the theoretical physics literature, see, e.g., [ |25| ] and references therein. Rigorous results are, however, 
rather scarce. An early example is [jig], where the authors get bounds on the free energy in the high 
temperature regime. Methods of stochastic analysis and large deviations are used in [ j34] ] to identify the 
limiting distribution of the partition function and also to obtain some information about the geometry of 
the Gibbs measure for small j3. More recent treatments of the high temperature regime using the very 



different methods are due to Talagrand [ |27| , see also [ |28| , Subsection 2.13]. The importance of the SK 
model with multidimensional spins for understanding the ultrametricity of the original SK model [ [26| l 
(which corresponds to d = 1 and ji being the Rademacher measure in the above notations) was emphasised 
in H). 

For the SK model, Guerra's scheme gave historically the first way to obtain the variational upper bound 
on the free energy in terms of the Parisi functional. The scheme is based on the comparison between 
two Gaussian processes: the first one being the original SK Hamiltonian ( |1.1[ ) and the other one being 
a carefully chosen GREM inspired process indexed by a G E^. The second important ingredient is a 
recursively defined non-linear comparison functional acting on the Gaussian processes indexed by a G E^y. 

The Aizenman-Sims-Starr (AS 2 ) scheme [jlj j2| gives an intrinsic way to obtain variational upper bounds 
on the free energy in the SK model. The scheme is also based on a comparison between two Gaussian pro- 
cesses. The first process is the sum of the original SK Hamiltonian X and a GREM-inspired process 
indexed by additional index space = N". The second one is another GREM-inspired process indexed by 
the extended configuration space x stf ' . The scheme uses a comparison functional defined on Gaussian 
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processes indexed by the extended configuration space equipped with the product measure between the 
original a priori measure and Ruelle's probability cascade (RPC) [^4j. The role of the comparison func- 
tional in the AS 2 scheme is played by a free energy functional acting on the Gaussian processes indexed 
by the extended configuration space. In p2[ ] Panchenko and Talagrand have reexpressed Guerra's scheme 
for the SK model using the RPC. 

Talagrand |30[ | using Guerra's scheme and the wealth of other ingenious analytical insights showed that 
the variational upper bound is also the lower bound for the free energy in the SK model. This established, 
hence, the remaining half of the Parisi formula. 

A particular case (d = 1, fl with bounded support) of the model we are considering here was treated 
by Panchenko in [pT|]. He used the techniques of [ p0[| to prove that in the case d = 1 upper and lower 
bounds on the free energy coincide (cf. (1.14) and (1.22) in this chapter). However, the results of [Ell 
Section 5 and the proofs of Theorems 2, 5 and 9] are based on relatively detailed differential properties of 
the optimal Lagrange multipliers in the saddle point optimisation problem of interest. These properties are 
harder to obt ain in m ultid imensional situations such as that we are dealing with here. In fact, as we show 
in Theorems 1.1 and 1.2, one can obtain the same saddle point variational principles without invoking the 



detailed properties of the optimal Lagrange multipliers. This is achieved using a quenched large deviations 
principle (LDP) of the Gartner-Ellis type. 

The most advanced recent study of spin-glass models with multidimensional spins was attempted by 
Panchenko and Talagrand in [|23||, where the multidimensional spherical spin-glass model was considered. 
The authors combined the techniques of [ f30| , [2l| ] to obtain partial results on the ultrametricity and also get 
some information on the local free energy for their model. 

Main results. In this paper, we prove upper and lower bounds on the free energy in the SK model with 
multidimensional spins in terms of variational inequalities involving the corresponding multidimensional 
generalisation of the Parisi functional (Theorems 1 1 . l| , [L2 , 5.1, 5.18). For this purpose, we generalise and 
unify the AS 2 and Guerra's schemes for the case of multidimensional spins, and employ a quenched LDP 
which may be of independent interest (Theorems [O and EO). Both schemes are formulated in a unifying 



framework based on the same comparison functional. The functional acts on Gaussian processes indexed 
by an extended configuration space as in the original AS 2 scheme. As a by-product, we provide also a 
short derivation of the remainder term in multidimensional Guerra's scheme (Theorem 5.4) using well- 
known properties of the RPC and the Bolthausen-Sznitman coalescent. This gives a clear meaning to the 
remainder in terms of averages with respect to a measure changed disorder. The change of measure is 
induced by a reweighting of the RPC using the exponentials of the GREM-inspired process]]. See [|2| for 
another approach in the case of the SK model (d = 1). 

We study the properties of the multidimensional Parisi functional by establishing a link between the 



functional and a certain class of non-linear partial differential equations (PDEs), see Propositions 6.1 
63] and Theorem p/2 . We extend the Parisi functional to a continuous functional on a compact space 



(Theorems 5.1, 5.2). We show that the class of PDEs corresponds to the Hamilton-Jacobi-Bellman (HJB) 



equations induced by a linear problem of diffusion control (Proposition 6.4 ). Motivated by a problem posed 
by [tnj], we show the strict convexity of the local Parisi functional in some cases (Theorem 5.4). 



We partially extend Talagrand's methodology of estimating the remainder term to the multidimensional 
setting (Theorem 5.4, Proposition |7.1[ Theorem 7.1). In the case of multidimensional Gaussian a priori 
distribution of spins we prove the validity of the Parisi formula (Theorem 1.3). 



We partially extend Talagrand's me thodology of estimating the remainder term to the multidimensional 
setting (Theorem 5A, Proposition [7. l| , Theorem 1_A). Though t he m ain technical problem of the method- 
ology in the general multidimensional setting remains (Remark ^5). In the case of the m ultid imensional 
Gaussian a priori distribution of spins we prove the validity of the Parisi formula (Theorem 1.3 ). 

Below we introduce the notations, assumptions and formulate our main results. The other results (men- 
tioned above) are formulated and proved in the subsequent sections. 

Assumption 1.1. Suppose that the configuration space E is bounded and such that G intconvE, where 
convE denotes the convex hull ofY.. 

The examples listed below verify this assumption: 



2 In d = 1 the latter fact was also known to the author of [[]], private communication. 
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(1) Multicomponent Ising spins. E = { — 1; l} d - the discrete hypercube. 

(2) Heisenberg spins. E={(7G[R' : ':||(t||2 = 1}- the unit Euclidean sphere. 

(3) E = {a G R d : ||ff||2 < 1} - the unit Euclidean ball. 

Remark 1.2. The boundedness assumption can be relaxed and replaced by concentration properties of the 
a priori measure. In Section |we will exemplify this in the case of a Gaussian a priori distribution. In 
general a subgaussian distribution will suffice. 

Consider the space of all symmetric matrices Sym(c/) = {A G R dxcl \ A = A*}. Denote 

Sym + (d) = {A G Sym(<f) | A h 0} , 

where the notation A >^ means that the matrix A is non-negative definite. We equip the space Sym(d) 
with the Frobenius (Hilbert-Schmidt) norm 

IMIf= L <v, MGSym(rf). 

M,V=1 

We shall also denote the corresponding (tracial) scalar product by (-,■}. For r > max{ || ct||| : a G E}, define 

W = {U G Sym(rf) \UhO, \\U\\ 2 < r} . 
We will call the set % the space of the admissible self-overlaps. In analogy to the usual overlap in the stan- 
dard SK model, we define, for two configurations, = (cr W ; o"i , . . . , <7^ ) G £#> ; = 1 , 2, the (mutual) 
overlap matrix Rn(g^\o^) G R rfxt/ whose entries are given by 

M^.^kv^^Lffffffff. ".veMnR (1.4) 

Fix an overlap matrix U G % . Given a subset "f C define the set of the /oca/ configurations, 

E w (r) = {c7GE w |^(a,cj)G y/}. 

Next, define the local free energy 

PN {Y) ee I log / e^^W). (1-5) 
W ^(r) 

We also define 

P (r) = p05,r)=jim^(n (i.6) 

where the existence of the limit follows from a result of Guerra and Toninelli [ |l9[ Theorem 1]. Consider a 
sequence of matrices =2 = {(QW g Sym(o')}^g such that 

= e (0) <Q [l) <...<Q {n+l) = U, (1.7) 

where the ordering is understood in the sense of the corresponding quadratic forms. Consider in addition a 
partition of the unit interval x = {x^Y]^, i.e., 

0=x <*\ < ■■■ <x n+ \ = 1. (1.8) 

Let {z^ }2 =0 be a sequence of independent Gaussian cZ-dimensional vectors with 

cov[ z «] = e (i+1) -e w . 

Given A G Sym(d), define 

X n+l {x,3,U,K) = \og f exp(V2j3(f;zit,cy> + (AtT,c7))dM(cy). (1.9) 
Define, for k G {«, . . . ,0}, by a descending recursion, 

Xt(x,^,U,A) = -logE zW [exp(xA +1 (x,^,f/,A))] (1.10) 

with 

X (x,^,E/,A) ee E z(0 ) [Xi(x,i2,E/,A)] , (1.11) 
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where E jk) [•] denotes the expectation with respect to the c-algebra generated by the random vector z®. 



Remark 1.3. Section 5.4 contains the more general framework of dealing with the recursive quantities 
1[ ) which in particular brings to light the links with certain non-linear parabolic PDEs. For these 
PDEs the recursion (1.2) is closely related to an iterative application of the well-known Hopf-Cole trans- 
formation, see, e.g., [Hy. 



Define the local Parisi functional as 

f{x,3,U ,A) = ~{A,U) - ^-j^ x k 

1 k=\ 



| G (*+1)||2_|| G «|,2 



2 +X (x,£,U,A). 



(1.12) 



Assumpt ion 1.2 (Hadamard squares). We shall say that a sequence, {Q^}" =1 , of matrices satisfies As- 
sumption 1.2 if 



( e <.») 



■■2 



-< 



-< 



(e w ) c %(e ( " +1) ) 



■■2 



(1.13) 



Remark 1.4. The above assumption on the matrix order parameters =2 is necessary only to employ the 
AS 2 scheme. In contrast, Guerra's scheme (Theorems 5.1 and 5.18) does not require the above assumption. 

One may verify that the matrices q and p in [|2i| T heorems 2.13.1 and 2.13.2] correspond to the matrices 
QW and Q<® of this paper (« = 1). (See also ( |1.25[ ) below.) Furthermore, a straightforward application of 
the Cauchy-Schwarz inequality shows that the matrices q and p actually satisfy Ass umption 1.2. We also 
note that in the simultaneous diagonalis atio n sc enar io in which the matrices in (L7) are diagonalisable in 
the same orthogonal basis (see Sections 63 and %2) this assumption is also satisfied. 

The first main result of the present paper uses the AS 2 scheme to establish the upper bound on the 
limiting free energy p(fi) in terms of the saddle point problem for the local Parisi functional ( 1.12 ). 



Theorem 1.1. For any closed set "f C Sym(<i), we have 



p(f)< sup inf f(x,3,A,U), 



(1.14) 



where the infimum runs over all x satisfying ( L8), all £2, satisfying both (L7) and Assumption L2 and all 
A G Sym(if ). 



We were not able to prove in general that the r.h.s. of (|1T4J) gives also the lower bound to the thermo- 
dynamic free energy. See, however, Theorem 1.3 for a positive example. 
To formulate the lower bound on (1.3) we need some additional definitions. 

Let the comparison index space be = N". Given a' 1 ', a' 2 ' £ £/, define 



,« (2) )) 



(1.15) 



(1.16) 



where q^(a^\ a' 2 ') is the (normalised) lexicographic overlap defined as follows 

, m m f°- a, (1 Va{ 2) 

ai (or ' a ( ') = 1 + < r i 

' [max\kell;n}n^:[a ( % = [aW]A, otherwise. 

Given a d x if -matrix M and p £ R, we denote by M 0p the d x ^-matrix with entries 

The matrix valued lexicographic overlap ( |1.15[ ) can be used to construct the multidimensional (d > 1 ) ver- 
sions of the GREM (see, e.g., [[|| and references therein for a review of the results on the one-dimensional 
case of the model). Here we shall need the following two GREM-inspired real-valued Gaussian processes: 
A = {A(a, Ct)} geZn ae&t an( l B = {B{a)}ae^ with covariance structures 



A(a 



)A(a^,a^) =2( J R((j( 1) ,(J^),G(« ll, ,a lzJ )>, 



E 



r(2)'i 
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Note that the process A can be represented in the following form: 

A{a,a) 



2 \ 1/2 n 

-) Z(M<*),°i), 



(1.17) 



where {A, = {A,(a)} cee ^}^ 1 are the i.i.d. (for different indices ;) Gaussian Revalued processes with the 
following covariance structure: for/S [1 ; JV] fl N, for all a^ l \ € £/ andallM,v€ [l;c/]nN assume that 
the following holds 

G(a (1) ,a (2) )„, v . 



Given t € [0; 1], we define the interpolating AS 2 Hamiltonian 

H,(a,a) = Vt{X(a)+B(a)) + VT~tA(a,a). (1.18) 
Next, we define the random probability measure 7Tjy G (E^ x &/) through 

where ^ = ^ (x) is the RPC [p4||. We denote by (cc)}a eiC / the enumeration of the atom locations of the 
RPC and consider the enumeration as a random measure on stf (independent of all other random variables 
around). Define the local AS 2 Gibbs measure ^v(f,x, i2,t/,'^ / ) by 

1 



f N (t,x,J2,U,r)\f] 



Z N (t,y) h N (f)x^ 



f{a,a)^ Ii ^^d% N {a,a) 1 



(1-19) 



where / : £yv x ^ —> R is an arbitrary measurable function for which the right-hand side of (1.19) is finite. 
For f Cf , define the AS 2 remainder term as 



9 N (f,x,£,U,r)®9 N (t,x,£,U,r) \\\R N {a w ,o {2) )-Q(a (l \a {2) )\\l 



dt. 



We define also the limiting AS remainder term 



&(x,£,U) = lim ]im& N (x,£ t B(U,e)) < 0, 



(1.20) 



(1.21) 



where B(U,e) is the ball with centre U and radius e. (The existence of the limiting remainder term is 
proved in Theorem 1.2.) 

The second main result of this paper uses the AS 2 scheme to establish a lower bound on (1.3) in terms 
of the same saddle point Parisi-type functional as in the upper bound which includes, however, the non- 
positive remainder term (1.21 ). In one-dimensional situations Talagrand [|o| and Panchenko [fTij], respec- 
tively, have shown that the corresponding error term vanishes on the optimiser of the Parisi functional. 



Theorem 1.2. For any open set V C Sym(c/), we have 



p(f)> sup inf \f{x,g,A,U)+@{x,£,U)], 



(1-22) 



where the infimum runs over all x satisfying (1.8), all A £ Sym(d), and all J3 satisfying both (1.7) and 
Assumption 1.2 

Remark 1.5. The comparison scheme of Guerra /[7^/ (see also more recent accounts /p2|/, / 7<S|/ and /p|/j 



is also applicable to our model and is covered by our quenched LDP approach, see Theorems 5.1 andp.l 



for the formal statements. Guerra's scheme seems to be more amenable (compared to the Aizenman-Sims- 
Starr one) for Talagrand 's remainder estimates /J?7J/, see Section ^ The scheme is based on the following 
interpolation 

H t {a,a) = VtX(a) + \f\~tA(a,a) (1.23) 

which induces the corresponding local Gibbs measure (1.19) and remainder term ( 1.20J ) by substituting 
( 1 . 1 8 |) with (1.23). Guerra's scheme does not include the process B and, hence, does not require Assump- 
tion L2. Recovering the terms correspondin g to <Pn(x, ty/)\B\ (see, ( 4.23 )) in the Parisi functional requires 
then a short additional calculation (Lemma 5.1). 
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Note that the results of Talagrand [ |28| , Theorems 2.13.2 and 2.13.3] imply that at least in the high 
temperature region (i.e., for small enough j3) the Parisi formula for the SK model with multidimensional 
spins is valid with n = 1 

p(P)=f(x,S*,0,U*)= sup inf f(x,4,A,U), (1.24) 
where the matrices g*' 2 ) = U* and Q*^' solve the following system of equations: 



d mf{x,£*,0,U*) = 0, M,ve [l;c/]nN, 

Qu,v 

d<i)f{x,B*,Q,U*) = 0, u,ve [l;rf]DlM. 



(1.25) 



Note that the system (1.25) coincides with the mean-field equations obtained in [ |28[ see (2.469) and 
(2.470)]. 

Let £ = R rf and fix some vector h G R d . Let fx G (£) be the finite measure with the following density 
(with respect to the Lebesgue measure A on £) 



/ \ 1/2 

d/i, , / detC \ / K . ., , 



(1.26) 



where C G Sym+(d). Note that, given m G R d and C G Sym + (rf) such that detC 7^ 0, the density ( |1.26[ ) with 
= Cm coincides (up to the constant factor exp (— l(Cm,m})) with the Gaussian density with covariance 
matrix C 1 and mean m. 

Remark 1.6. It turns out that only matrices C with sufficiently large eigenvalues will result infinite global 
free energy, cf. Lemma 8.8. The local free energy is, in contrast, always finite, see Lemma 8.7 and Theo- 
rem |/.4 

Consider the function / : (0 : +°°) 2 



f(c,u) = 



fj3V 



R given by 
log cm — cu + 1 , 



ug(0;$] 



\(2V2j3- C ) M + log^-I(l+log2) 



«e(#;- 



(1.27) 



The following result shows that, at least, in the highly symmetric situation (1.26) with h = the multidi- 
mensional Parisi formula indeed holds true (see Lemma 8.7 for an explanation why the result is indeed a 
Parisi formula). 



Theorem 1.3. Let jX satisfy (1.26) with h = 0. Assume that the matrices U and C are simultaneously 
diagonalisable in the same basis. Denote by {c v G R+}y = i and {u v G R+}y = j the eigenvalues of the 
matrices C and U, respectively. Moreover, assume that min,, u v > and min v c v > 0. 
Then we have 

d 

lim lim p N (L N (B(U ,e))) = Yf(c v ,u r ). 

e|+0WT+~ ,fl 

Remark 1.7. Close results have previously been obtained in the case of the spherical model in /]23|/, from 
where we borrow the general methodology of the proof of the Theorem 1.3 As noted in /]2i[/, another 
more straightforward way to obtain the Theorem 1.3 is to diagonalise the interaction matrix G and use the 
properties of the corresponding random matrix ensemble. 

Organisation of the paper. The rest of the present paper is organised as follows. In Section ||] we record 
some basic properties of the covariance structure of the process X and establish the relevant concentration 
of measure results. The section contains also the tools allowing to compare and interpolate between the 
free energy-like functionals of different Gaussian processes. In Section pi we derive a quenched LDP of 
the Gartner-Ellis type under measure concentration assumptions. Section^ contains the derivation (based 
on the AS 2 scheme) of the upper and lower bounds on the free energy of the SK model with multidimen- 
sional spins in terms of the saddle point of the Parisi-like functional. In Section || we employ the ideas 
of Guerra's comparison scheme in order to obtain the upper and lower bounds on the free energy and we 
also get a useful analytic representation of the remainder term. In Section || we study the properties of 
the multidimensional Parisi functional. Section [7] contains the partial extension of Talagrand's remainder 
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term estimates to the case of multidimensional spins. In Section || a case of Gaussian a priori distribution 
of spins is considered and the corresponding local Parisi formula is proved. In the appendix we prove the 
almost super-additivity of the local free energy, as an application of the Gaussian comparison results of 



Subsection 2.3 



2. Some preliminary results 



2.1. Covariance structure. Our definition of the overlap matrix in (1.4) is motivated by the fact that, as 
can be seen from a straightforward computation 

e[x n (g^)X n (g^)] = £ (R N (o( 1 \o^ v y=\\R N (oW,oW)\\j, (2.1) 

if,V=l 

that is, the the covariance structure of the process Xn(g) is given by the square of the Frobenius (Hilbert- 
Schmidt) norm of the matrix Rn{g^ 1 \ G^). The basic properties of the overlap matrix are summarised in 
the following proposition. 



Proposition 2.1. We have, for all G^ , G^ , G € E^v, 

(1) Matrix representation. R n (g { - 1 \ <7 (2) ) = jj ( cj(1) ) 

(2) Symmetry #1 . R u / (a< l \ G^ ) = (g^ , a< 1} ). 



r(2) 



(3) Symmetry #2. R"/(g, g) = R^"(g, g). 

(4) Non-negative definiteness #1. Rn(g, g) >: 0. 

(5) Non-negative definiteness #2. 

~R n (gW,gW) Rjf(aW,a&>) 
_R n (gW,gW)* R n {gW,gW) 



^0. 



(6) Suppose U = R n (g^ ,G^) =R n (g ( - 2 \g {2 1), then 

||*(a(V (2) )lll<l|C/|||. 

Proof. The proof is straightforward. □ 

2.2. Concentration of measure. The following concentration of measure result for the free energy is 
standard. 

Proposition 2.2. Let (E, &) be a Polish space. Suppose jj, is a random finite measure on E. Suppose, 
moreover, thatX(G), G G E is the family of Gaussian random variables independent of pi which possesses 
a bounded covariance, i.e., 

there exists K > such that sup \Cov(X(g {1) ),X(g {2) ))\ < K. (2.2) 

<tU),ct(2) g i 

Assume that 

f(X)=l0gj^ X ^dfl(G)<oo. 

Then 

t 2 \ 



F{\f(X) - E[f(X)}\ > t} < 2exp (- — J 

Remark 2.1. An analogous result was given in a somewhat more specialised case in /|2~/|/. 

Proof. This is an adaptation of the proof of [|28[ Theorem 2.2.4] . We can not apply the comparison Theorem 



2.5 directly, so we resort to the basic interpolation argument as stated in Proposition 2.1. For j = 1,2, let 
the processes Xj(-) be the two independent copies of the process X(-). For t 6 [0; 1], let 

Xjj = y/iXj + \f\~tX 



and 



Fj(t) = \og j exp(Xj,(G))d^G). 
Jn 
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For s € R, let 

Hence, differentiation gives 



(p s (t) = E[exp (s(F l -F 2 ))}. 



(p s (t) = sE [exp 0(F] - F 2 )) (Fi - Ft)] 
(the dots indicate the derivatives with respect to t ) and also 



(2.3) 



J (t- l l 2 Xj{a) - (1 -t)- l l 2 X{a)) e X p(X 7V ((T))dM(cT). (2.4) 



Now, we substitute (2.4) back to (2.3) and apply Corollary 2.1 to the result. After some tedious but ele- 
mentary calculations we get 



<p s (t) =s 2 E 



exp(i(Fi-F 2 )) (J^expX h ,(a)dfi(a)J^expX ltt (a)d^{a] 



^Cov(X((7 (1 )),X((J (2 )))exp(x L ,((7( 1) )+X2, ; ((7 (2) ))d^((7( 1) )d^(c7( 



Thus, thanks to (2.2), we obtain 

(j) s (t)<Ks 2 (p s (t). 

The conclusion of the theorem follows now exactly as in the proof of [M Theorem 2.2.4]. 



□ 



We now apply this general result to the our model and also to the free energy-like functional of the 
GREM-inspired process A. 

Proposition 2.3. Suppose E C B(0, r), for r > 0. For Q. C T.n, denote 

4 K (P,Cl) ee log^exp (VNpX N ((T)) d^ N (a), 



and 



P^(j3,0)EElogj[ ix ^exp(i3V2£(A,(a) J cr i ))d^(tT ! a). 



Then, for all £1 C £/v, we have 

(1) For any t > 0, 

P{|^,fl) - E [^(J3,a)] | > r } < 2 exp 

(2) For any t > 0, 

P{|/^ £M (/3,n) - E [P™(j3,i2)] | > t) < 2exp 



8/3 ViV 



(2.5) 



(2.6) 



Proof. (1) We would like to use Proposition [2.2[ By ( |2.1[ ) and the Cauchy-Bouniakovsky-Schwarz 
inequality, we have, for all WeN, ffW,ff™ £ that 

Cov(^(a( 1 ),a( 2 ))) = ||^(cT( 1 ),cT( 2 ))|| 2 = -l £ (a/ 1) ,a] 1) )<a i (2) ,cTf ) (2.7) 



>J=i 



Hence, for all N € N and all subsets £2 of E#, we obtain 

sup \Cow(X(a il) ),X(a {2) ))\<r 4 . 



Thus (2.5) is proved. 



Anton Bovier and Anton Klimovsky 
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(2) We fix an arbitrary N G N, 0^,0® G L N , a {l) ,a^ G ^. We have 

Cov(A(a (1) ,a (1) ),A(cj( 2 ),a( 2 ')) = E[A((j (1) ,a (1) )A(a (2) ,a (2) ) 

= £<fi(«(V (2) M ( V/ 2) >. 



Bound (2.7) implies that, for any U G ^, we have ||C/j| 9 < r 2 . Since g(a (1 ',a' 2 ') G we obtain 



|(G(a« I aP))o / (1 U (2 ))|<||e(a(^a^0ll2||cTr ; ll2|^ril2 
<||G(a« a (2) )lk 2 <r 4 . 



Therefore, using Proposition |2.2], we obtain (2.6). 



□ 



2.3. Gaussian comparison inequalities for free energy-like functionals. We begin by recalling well- 
known integration by parts formula which is the source of many comparison results for functionals of 
Gaussian processes. 

Let F : X — ► IR be a functional on a linear space X. Given x EX and e E X, a directional (Gateaux) 
derivative of F at x along the direction e is 



d x ^ e F(x) = d t F(x + te) . (2.8) 

r=0 

With this notation the following lemma holds. 

Lemma 2.1. Let {g(i)}iei be a real-valued Gaussian process ( the set I is an arbitrary index set), and h be 
some Gaussian random variable. Define the vector e G IR as e (/) = E [hg{i)\, i El. LetF-.R'^R such 
that, for all f E R 1 , the function 

R3t^F(f + te)ER (2.9) 

is either locally absolute continuous or everywhere differentiable on R. Moreover, assume that the random 
variables hF{g) and d g ^ e F{g) are in L 1 . 
Then 

E[hF(g)} = E[d g ^ e F(g)}. (2.10) 

The previous proposition coincides with [ fTl} Lemma 4] (modulo the differentiability condition on (2.9) 
and the integrability assumptions which are needed, e.g., for [|[ Theorem 5. 1.2]). 

The following proposition connects the computation of the derivative of the free energy with respect to 
the parameter that linearly occurs in the Hamiltonian with a certain Gibbs average for a replicated system. 

Proposition 2.4. Consider a Polish measure space (£,6) and a random measure fi on it. Let X = 
{X((7)} CTG £ and Y = {Y((j)} ae z be two independent Gaussian real-valued processes. For u E R, we 
define 

H u (a) = uX(a) + Y(a). 
Assume that, for all u E [a,b] <s IR, we have 

J exp(//„(c7))djLi((7) <°o,J X(a)exTp(H u (a))dfi(a) <°° 

almost surely, and also that 



log / exp(H H (cr))dju(c7) 



Then we have 



du 



E 



log 



d M (CT) 



: uE p{u) ® &(u) [VarX(a) - E [X(a),X{x)}}} , 



where (m) is the random element of \M\ (£) which, for any measurable f : L — * R , satisfies 

1 



Z(u) 



/(cr)exp( J f/ M (cj))d J L i (cT). 



12 



The Aizenman-Sims-Starr and Guerra's schemes for the SK model with multidimensional spins 



Proof. We write 



du 



-flog / e H "(°>dii(<J) = I X(a)^-d^(a), 



Ma) 



UP) 



(2.11) 



where Z„(J3) = / e^^d^a). The main ingredient of the proof is the Gaussian integration by parts 
formula. Denote, for T€l, e(z) = E [X{a)H u (x)}. By fl2.10|), we have 



X(a) 



e H u (a) 



= E 



dx 



Ma) 



/e^(*)d/i(T) 



(2.12) 



Due to the independence, we have 



E[Z(ct)H h (t)] = mE[Z(c7),Z(t)]. 



Henceforth, the computation of the directional derivative in (2.12) amounts to 



d_ 

dt 



Jl u {a)+tuN?a{a) 



| e H„(T)+f«Cov(cr,T) ( j^( T ) 



= ( y e // "( fj) d J u((j)J ^«VarX((7)e // "( tJ ' y e^d^r) 
-e H " (<7 > / uCow[X{o),X{x)]e. H ' M dil{%) 



Substituting the r.h.s. of (2.13) into (2.1 1 ), we obtain the assertion of the proposition. 



(2.13) 

□ 



The following proposition gives a short differentiation formula, which is useful in getting comparison 
results between the (free energy-like) functionals of Gaussian processes. 

Proposition 2.5. Let {X{c)) ae z> oeL be two independent Gaussian processes as before. Set 

H,(a) = VtX(a) + VT~^tY(a). 

Assume that 

J e H ' {a) d^i{a) <o°,j X(a)e Hl{a) dii(a) < °°, 
J Y{a)e H ' ia) dii{a) < «, 

almost surely, and also that, for all t G [0;1], 



log / e M > [a >dn(<7) 



< oo. 



Then we have 



log J e x{a) d^i((j) = E log J e Y ^dii(a) 

1 [ W(0<8)Sf(/) [fvarX(a (1) )- VarF(ffW) 

2 Jo LV 



Cov 



X(c7 (1 )),X(a (2) ) -Cov F(CT«),y(<7 



df, 



where §f (?) is f/ie random element o/^#i(E) which, for all measurable f : E — > IR, satisfies 

1 



^(0[/] 



Z(f) A 



/(CT)«p(fli(ff))dM(ff). 



(2.14) 



(2.15) 



Proof. Let us introduce the process 



Hence, 



W M , v (cr)EE M X(a)+v7((T) 



(2.16) 
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Thus 



log 



Ma 



d M (C7) 



]_d_ 
2 \ ~\/t ~du 
1 



log 



E 



>d M (c7) 



log / e w "^ a 'dii(a) 



Applying Proposition ^4 and / df to the previous formula, we conclude the proof. 



□ 



3. Quenched Gartner-Ellis type LDP 



In this section, we derive a quenched LDP under measure concentration assumptions. Theorems 3.1 



and 3.2 give the corresponding LDP upper and lower bounds, respectively. The proofs of the LDP bounds 
will be adapted to get the proofs of the upper and lower bounds on the free energy of the SK model with 
multidimensional spins. However, they may be of independent interest. 

Note that the existing "level-2" quenched large deviation results of Comets JlO| ] are applicable only 
to a certain class of mean-field random Hamiltonians which are required to be "macroscopic" functionals 
of the joint empirical distribution of the random variables representing the disorder and the independent 
spin variables. The SK Hamiltonian can not be represented in such form, since the interaction matrix 
consists of i.i.d. random variables. Moreover, it is assumed in [|o[ that the Hamiltonian has the form 
Hn(°) = NV(o), where {V(a)} a eZ N is a random process taking values in some fixed bounded subset of 
IR. Since the Hamiltonian of our model is a Gaussian process, this assumption is also not satisfied, due to 
the unboundedness of the Gaussian distribution. 

3.1. Quenched LDP upper bound. The following assumption will be satisfied for the applications we 
have in mind. As is clear from what follows, much weaker concentration functions are also allowed. 

Assumption 3.1. Suppose {Qn}n=i ' s a sequence of random measures on a Polish space , X). Assume 
that there exists some L > such that for any Q^-measurable set A C X we have 



P{|logftv(A) - E[logftv(A)]| > t} < exp 



t 

In J 



(3.1) 



Note that Assumption 3.1 will hold in the cases we are interested in due to Proposition 2.2 



Lemma 3.1. Suppose {Qn}n=\ i s a sequence of random measures on a Polish space (JT,X) and for 
{A r C X : r G {1,...,/?}} is a sequence of Qn -measurable sets such that, for some absolute constant 
L > and some concentration function TjN{t) : IR+ — > R+ with the property 



rj N (tN)dt 



0. 



we have 



Then we have 



P{|log<av(Ar)-E[log&v(A r )]|>?} < IJaKO- 



lim -E 



logftv(L> 

r= 1 



- max E[logQ N (A r 
re{l,..., P } 



= 0. 



(3.2) 



(3.3) 



(3.4) 



Remark 3.1. As is easy to extract from Assumption 3.1 we will apply this result in the very pleasant 
situation, where 

,2 



YN{t) 



However, our subsequent results hold for substantially worse concentration functions satisfying (3.2). 



Proof of Lemma 3.1. First, (3.3) gives 
P 



max |log<2jv(Ar)-E[log(MA r ) 
re{l,...,p} 



>t><pr} N (t). 
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Since, for a,b G R p , the following elementary inequality holds 



maxfl,- — maxb r 



< max|a r — b r \, 



we get 



max \ogQ N (A r )- max £[logQ N (A r )} 

re{l,...,p} re{l,...,p} 



>t} < P rj N (tN). 



The last equation in turn implies that 



max loggAr(A r )- max E[logQ N (A r )] 

re{l,...,p} re{l....,p} 



<p ( 7] N {tN)dt, 
Jo 



and the r.h.s. of the previous formula vanishes as N | +°° due to (p.2|). 

Let Qn G {3s), N G N be a family of random measures on («^T,X). Define the Laplace transform 

L*(A) = / e^ A >d&v(;c). 



(3.5) 

□ 



Suppose that, for all A G R , we have 



/(A) = lim— E[logLjv(A)l G R = R U {-■», +00} . 



Define the Legendre transform 



Define, for 5 > 0, 



Lemma 3.2. Suppose 



Then 



r(*)=inf[-<x,A)+/(A)]. 

A 



/|(x)=max<jr(x) + 5,-- 



G int^(f) = int{A : /(A) < +00} . 



(3.6) 



(3.7) 



(3.8) 



(3.9) 



(1) The mapping I* (•) iJ-tRis upper semi-continuous and concave. 

(2) For a// M > 0, 

{.1 £ 1" : /* (jc) < M} is a compact. 
Proof. (1) Since, for all A G the linear mappings 

x^-(A,x)+I(A) 

are obviously concave, the infimum of this family is upper semi-continuous and concave. 
(2) See, e.g., [||] for the proof. 

Theorem 3.1. Suppose that 



□ 



(1) The family {Qn} satisfies condition (3.4). 

(2) Condition (3.6) is satisfied. 

(3) Condition (3.9) is satisfied. 

Then, for any closed set C R d , we have 



lim-E[log&v(r)]<su P r(x). 



(3.10) 
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Proof. (1) Suppose at first that V is a compact. 

Thanks to (3.7), for any x G St ' , there exists A(x) <G St such that 

-<x,A(x))+/(A(x))</|(x). (3.11) 

For any x G Sf, there exists a neighbourhood A (x) C St of x such that 

sup (y — x,A(x)) < 8. 

y eA(x) 

By compactness, the covering \J xe ^A(x) D "V has the finite subcovering, say Ur=i^( x r) 3 
Hence, 



1 

at 



iogftv(r)<-io g ( U&(^W) 



v.r=l 



Applying condition (3.4), we get 



lim — E 



max log (A (x,-)) — max E 

re{l,...,p} rG{l,...,p} 



2V 



logftv(A(* r )) 



<0. 



(3.12) 



(3.13) 



By the Chebyshev inequality, 

Qn(A(x)) <Q N {y£S£: (y-x,A(x)) < 8} 
<e- SN J e N< J-*M*)) d Q N (y) 
^ e -SN e -N( x ,M,)) LN{A{x) y 



(3.14) 



Hence, (3.14) together with(3.11)yields 
km -J-E[log&v(A(x r ))] < lim 



■(x r ,A(x r )) + -logIiv(A(x r )) 



= -(x r ,A(x r ))+/(A(x r ))-5 
<4(x,)-5. 



(3.15) 



Combining fl3.12| ), ( p.l3[ ), ( p.l5[ ), we obtain 

1 



lim - E [log Q N (Y)]< max £ (x r ) - 5 

N]+oolS re{l,...,p] 

< sup/g(x) — 5. 

Taking 5 J. +0 limit, we get the assertion of the theorem. 
(2) Let us allow now the set Y to be unbounded. We first prove that the family Qn is quenched 
exponentially tight. For that purpose, let 

R N (M) ee Ie [log<av( Jf \ [-M;M] rf ) 



and denote 



We want to prove that 



R(M) ee lim Rn(M). 



lim Z?(M) 



(3.16) 



Fix some u G Suppose 5„ p G {0, 1} is the standard Kronecker symbol. Let e u G R d be 

an element of the standard basis of R d , i.e., for all p G {1, ... ,d}, we have 

Thanks to the Chebyshev inequality, we have 



giv{x H < — M} < e 



-NM 



-N(x,e u ) 



dQ N (x), a.s. 



(3.17) 
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Now, we get 



Ln(A c — e u ) 



Hence, combining (3.17) and (3.18), we obtain 
1 



^E[log&v{x„ < —M}] < -M+I N {A e -e u )-l N {K)- 



Using the same argument, we also get 
1 



N 



E[logQ N {x u >M}] < -M + I N {A e +e u ) -I N (A e ). 



We obviously have 

1 r / d 
Rn(M) < -E [logftv ( (J {{x u < —M} U {x u > M}) 

€ w=l 
4h to (3.21 ), we get 

EE i E [log Q N ( |J ({x„ < —M] U {x„ > M}) ) 

{E[log&v({x„ < -Jlf})],Epog(2v({x, > Af})]} 



(3.18) 



(3.19) 



(3.20) 



(3.21) 



max max-; 

«G{1 d} 



<0. 



(3.22) 



Applying ( |3.19| ) and ( |3.20| ) in ( |3.22[ ), we get 



lim — E 

N]+o°N 



log ( U ({*« ^ - M > u i A » ^ M » , 

u=\ 

<—M — I(A e )+ max max{/(A e — e„),/(A e + e„)} . 
«e{l,...,d} 



(3.23) 



The bound (3.23) assures (3.16). Now, since we have (with the help of (3.4) and (3.10)) 
Gin" -J-E[log&v0n]< K ^-E[loge iv ((rnhM;Mf)U(jr\[-M;M] d )) 



<max{ sup I*{x),R{M) 



(3.24) 



the assertion of the theorem follows from ( 3.16 ) by taking the lim^+oo in the bound ( 3.24 ) 

3.2. Quenched LDP lower bound. Suppose that, for some A£R d and all N G N, we have 

J e N ^dQ N (y) < +oo. 

Let Qn.a G ^(3£) be the random measure defined by 

Qn.a(A)= [ e N ^dQ N (y), 

J A 

for any Qn measurable A C 2£ '. 



□ 



(3.25) 



Lemma 3.3. Suppose the family of random measures Qn satisfies the following assumptions. 

(1) Measure concentration. For all N G H there exists some L > anc/ rjw : R+ — > R+ smc/z f/zaf, for 
any Qn -measurable set A C we have 

F{\\ogQ N (A)-E[logQ N (A)]\>t}<r lN (t). 

Assume, in addition, that, for some p > 0, the concentration function satisfies 

N p f r] N (Nt)dt >0. (3.26) 

Jo N1+<» 
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(2) Tails decay condition. Let 



C(M) = {xe$? :\\x\\ <M}. 



There exists p <E N such that 



lim lim / + v\-\ogQ NA (&\C(N p ))>-K+t\dl = 0. 

K]+o°N]«.J() [N ' 



(3.27) 



(3) Non-degeneracy. The family of the sets {B j C 3£ : j G {1,...,^}} satisfies the following condition 



there exists some jo e {1, . . . ,q} such that lim — E \ogQfj A (B j ) 



N 



(3.28) 



77!e«, /or a«y A £ R , we have 



1 



lim — E 



log 



G*,a(U b /! 



;=1 









) — max E 


log&v,A0B/) 




> Ml,-,?} 





<0. 



(3.29) 



Remark 3.2. 77ie polynomial growth choice of M = A^v = made in assumptions (3.27) anc/ (3.26) is 
made for specificity. Inspecting the following proof one can easily restate the conditions (3.27) and (3.26) 
for general dependencies. Effectively, the growth rate ofMpj is related to the covering dimension of the 
Polish space (St^,X). 



Proof of Lemma [3. 4 We fix some j G {l,...,g}. Take an arbitrary e > 0, M > and denote Jm,e 
Zn [— ||A||M/e; ||A[|M/e]. Consider, for i € 7m,e, the following closed sets 



We get 



We have 



Denote 



and 



Ajj = {xe Bj : (j - l)e < (A,jc) < je}. 



^ log Qn.a ( U Bj) < i log Q N . A ( ( [J Bj n C(M)) U ( \ C(M) ) 
•/ =1 , . > =1 



< 



- max {^max logQ NA (Bj n C(Af)), 

ioge,, A (^\J(M))}^ log(9+1) 



2V 



^logG W ,A(B,nC(M)) < I log ( £ c m Q N (A itJ )) 



< max 

i£{l,..,p} 



ie + ^logQ N (A L j) 



log(cardy M ,e) 
iV 



ajv(g)= max max /£ + — :loggw(A/ ,) 

y'6{l,..,?}i'€/M,£ V N 



ji N = max E 

;£{!,..,?} 



loge^A^i) 



(3.30) 



(3.31) 



j3jv(e) = max E 
Mi.-.?} 



max i£ + -logQ N (Aij) 

ieJ M .F \ Jv 



y N (M) = -logQ NA (^\C(M)). 



18 



The Aizenman-Sims-Starr and Guerra's schemes for the SK model with multidimensional spins 



We also have 



^logQ NA (Bj) > ^logQ NA (Bjr\C(M)) 



> max 



max 



i£ + -\og Q N {Aij] 



Due to condition ([[]), we have 

p{ \a N (e)-p N (e) >f| < rj N (tN)qcaidJ M , £ . 
We put M = Mn = N p , and we get 

card7 M ,£ <2||A||M/e + l 
<2||A||JV p /£ + l. 

Let 

X N {M, e) = max{7 W (M), a N {e)} - fi N , 

then we have 

P{X N (K,e) >t}< P{yn(M) >p N + t} + P{a N {e) > p N + t}. 
Due to property ( 3.28| ), there exists K > such that we have 

P{Y N (M) >p N + t}< P{y N (M) >-K+t}. 



Thanks to ( 3.32 ), we have 



P{a N (e) >p N + t}< P{a N (e) > p N ( e )+t- e}. 



For t > e, we apply ( |3.33| ) and ( |3.34| ) to ( |3.37[ ) to obtain 

P{a N (e) >p N + t} < (2\\A\\N p /e + l)qri N (tN). 



Combining ( 3.3C ) and ( 3.31 ), we get 



E 



l0g&V,A (U S ./ 

L ;'=i 
<E[X N (M,e)} + 

Now, and (gj|) imply 



max E 



log&v, A (B;) 
log( 9 +l) log(2||A||JV"/e + l) 



EKiv(M,e)] < / P{Xjv(M,e) >t}dt 
Jo 

< J P{X N (M,e) > t}dt + e 
+ (2\\A\\N"/e + l)q ( + T] N (tN)dt + e. 



(3.32) 



(3.33) 



(3.34) 



(3.35) 
(3.36) 
(3.37) 
(3.38) 



(3.39) 



(3.40) 



Therefore, taking sequentially \im^ +0! „ lim^| +0O and lim £ | + o in ( 3.40| ), and using (3.26), we arrive at 



lim£\X N (M,e)} < 0. 



Bound (3.41) together with (3.39) implies the assertion of the lemma. 



(3.41) 



□ 
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Let Qn.a be the (random) probability measure defined by 

Lemma 3.4. S uppose that the measure Qn satisfies the assumptions of the previous lemma. 
Then ( 3.29 ) is valid also for Qn.k- 

Proof. Similar to the one of the previous lemma. □ 

Remark 3.3. Recall that a point x G 5£ is called an exposed point of the concave mapping I* if there exists 
A G R d such that, for ally G SC\ {x}, we have 

r~(y)-I*(x)<(y-x,A). (3.42) 

Theorem 3.2. Suppose 



(1) The family {Qn ■ N G N} C satisfies the assumptions of Lemma 3.3 

(2) §f C 3£ is an open set. 

(3) 7^ C S>{I*) is the set of the exposed points of the mapping I*. 

(4) Condition (3.9) is satisfied. 



Then 



M -E[log&v(Sfn^)]>su P r(x). 



N] + 

Proof. Let B(x, e) be a ball of radius e > around some arbitrary x G SC . It suffices to prove that 

1 



Indeed, since we have 



lim lira - E [log ftv (B (x, e ))]> I* (x) . 



Qn{^)>Qn(B{x, £ )), 



(3.43) 



(3.44) 



(3.45) 



apply ing ^ log (-), ta king the expectation, taking firrw +00 , e J, +0 and taking the supremum over x G Sf in 
( p74^ ), we get ( P~43b . 

Take any .tg^nl. Then we can find the corresponding vector A e = A e (x) G R orthogonal to the 
exposing hyperplane at the point x, as in (3.42 ). Define the new ("tilted") random probability measure Qn 
on R d by demanding that 



00 = 



1 



JV<y,A,) 



dg w L N (A e ) 



(3.46) 



Moreover, we have 



iE[log&vOB(*,£))] = ^E 



log / dQ N (y) 

JB(x,e) 

^E[logZw(A e )] + iE 



B(x,e) 



dftv(y) 



Hence, 



> ^E[\ogL N (A e )} - (x,A e ) - e\\A e \\ 2 + ^E[\ogQ N (B(x,e))] 



lim ljml-E[logQ N (B(x,e))] > [-(x,A e ) +I(A e )} + lim Mv r E [log ftv (#(*,£))] 



Since we have 



■{x,A e )+I{A e )>I*{x), 



in order to show (3.44) it remains to prove that 

1 



lim lim -E [log Q N (B(x,e))] = 0. 
el+OATfooA' 



(3.47) 
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The Laplace transform of Qn is 



Hence, we arrive at 



Moreover, we have 



L N {A) = 



L N (A + A e ) 



L N {A e ) 
/(A)=/(A + A e )-/(A e ). 

!*(x)=r(x) + {x,A e )-I(K)- 



(3.48) 



By the assumptions of the theorem, the family Qn satisfies the assumptions of Lemma 3.3 Hence, due to 
Lemma 3.4, the family Qn satisfies (3.4). Thus we can apply Theorem 3.1 to obtain 



lim — E 

N]+°° N 



\ogQ N (R d \B(U,e)) < sup !*(y). 

J ye<W\B{x,E) 



Lemma implies that there exists some xq € S£ \ B(x, e) (note that xq ^ x) such that 

sup f (y)=t*(x ). 

yeStr\B(x,e) 



Since A e is an exposing hyperplane, using (3.48), we get 

!* (x )=I*(x ) + {xo,Ae}-I(K) 

< [I*(x Q ) + (x ih Ae)] - [/*(*) + (x,K)} < 0, 
and hence, combining ( 3.49[ ) and (3.50), we get 



lim — E 



logQ N (R d \B(x,e)) 



<0. 



Therefore, due to the concentration of measure, we have almost surely 

Urn" ±-logQ N (R d \B{x,e))<0 

N]+<*>N 

which implies that, for all £ > 0, we have almost surely 

lim Q N (R d \B(x,e))=Q, 



and (3.47) follows by yet another application of the concentration of measure. 



(3.49) 



(3.50) 



□ 



Corollary 3.1. Suppose that in addition to the assumptions of previous Theorem 3.2 we have 

(1) /(■) is differentiate on int^(7). 

(2) Either 9(1) = % or 



Then £"(!*) = R d , consequently 



lim ||W(A)|| = +oo. 
A->dJ>(/) 



lim iE[loge^(^)]>su P r(x). 



Proof. The proof is the same as in the classical Gartner-Ellis theorem (see, e.g., [D_3J) 



□ 
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4. The Aizenman-Sims-Starr comparison scheme 



In this section, we shall extend the AS scheme to the case of the SK model with multidimensional 
spins and prove Theorems [Q] and [Eg, as stated in the introduction. We use the Gaussian comparison 
results of Section 2.3 in the spirit of AS 2 scheme in order to relate the free energy of the SK model with 
multidimensional spins with the free energy of a certain GREM-inspired model. Comparing to [[!]], due to 
more intricate nature of spin configuration space, some new effects occur. In particular, the remainder term 
of the Gaussian comparison non-trivially depends on the variances and covariances of the Hamiltonians 
under comparison. To deal with this obstacle, we use the idea of localisation to the configurations having a 
given overlap (cf. ( |1.5[ )). This idea is formalised by adapting the proofs of the quenched Gartner-Ellis type 
LDP obtained in Section 

4. 1 . Naive comparison scheme. We start by recalling the basic principles of the AS 2 comparison scheme 
(see, e.g., Chapter 11]). It is a simple idea to get the comparison inequalities by adding some additional 
structure into the model. However, the way the additional structure is attached to the model might be 
suggested by the model itself. Later on we shall encounter a real-world use of this trick. Let (£,G) and 
be Polish spaces equipped with measures ji and respectively. Furthermore, let 

X = {X(a)} aeZ ,A = {A(a, a)} aeL . ,B = {B{o)} aea ? 
be independent real-valued Gaussian processes. Define the comparison functional 

4>[C] = Elog/ e cia < a) d(n®£){o,a) , (4.1) 

where C = {C(a, a)} a ez is a suitable real-valued Gaussian process. Theorem 4.1 of [|[] is easily under- 

stood as an example of the following observation. Suppose <P[X] is somehow hard to compute directly, but 
<t>[A] and <P[B] are manageable. We always have the following additivity property 



&[X+B] =$[X]+&[B}. 

Assume now that 

<P[X + B] < 4>[A] 

which we can obtain, e.g., from Proposition |2.5|. Combining (4.2) and (4.3), we get the bound 

<& [X] < 4> [A] - 4> [B] . 
4.2. Free energy upper bound. Let "V c Sym(rf) be an arbitrary Borell set. 
Remark 4.1. Note that is closed and convex. 
Let 

£ W (T) = {a e Z N 
= {a G E/v 

Let us define the local comparison functional < t , Ar(x, as follows (cf. (4.1)) 

1 



(4.2) 



(4.3) 



(4.4) 



R N (<?,<J) € fnf}. 



<f w (x,r)[c] 



-E 



N 



\ogn N 



(4.5) 



(4.6) 



where C = {C(<J, a)} a ez is a suitable Gaussian process. Let us consider the following family (Af £ N) of 
random measures on the Borell subsets of Sym(c/) generated by the SK Hamiltonian, 



P N (T) = [ c^^d^ia), 



and consider also the following family of the random measures generated by the Hamiltonian A(a, a) 

N 



exp J3 VnY (Ai(a), <r f ) dn N {a, a), (4.7) 

where the parameters £2 and U are taken fro m the definition of the process A(ot) (cf. (L7)). The vector x 
defines the random measure t, € (s>/) (cf. (L8)), and, hence, also the measure Kn G ^#(E x &/). 
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Remark 4.2. To lighten the notation, most of the time we shall not indicate explicitly the dependence of 
the following quantities on the parameters x, i2, U. 



Consider (if it exists) the Laplace transform of the measure (4.7) 

L N (A)^ [ c N ^6P N (U). 

J"?/ 

Let (if it exists) 



/(A) = lim-E logZw(A) 

iVT~ N L 



Define the following Legendre transform 



I*{U)= inf -<{/,A)-4>A,(x,r)[B]+/(A) 

JCG^'(l.l), L 
AGSym(d) 



(4.8) 



(4.9) 



(4.10) 



Denote, for 8 > 0, 



Let 



/|(I7)=max<jr(C/) + 5,-- 



P {-r)= iim -E[io g /v(r)]. 



(4.11) 



Remark 4.3. Note that the result o//^/ assures the existence of the limit in the previous formula. 

Lemma 4.1. We have 

( 1 ) The Laplace transform ( X. 8 ) exists. Moreover, for any A £ Sym (of), we have 

e W)dP w ([/) 

exp (n(A,R n (o,o)) +pVNX{(jj) dn m (a), 



(4.12) 



= / exp(Ar(A,i? iV (c7,C7))+j3\/iV£(A i (a),C7 i ))d^ v (tT,a). (4.13) 

(2) 77ie quenched cumulant generating function (4.9) erisfs in f/ie | °° /wiit /or any A £ Sym(af). 
Moreover, for all N £ N, we /lave 



7 W (A) = ~EpogLtf(A)] =X (^^,A,C7), 



(4.14) 



f/zaf /i In(') in fact does not depend on N. 



Proof. (1) We prove ( 4.13 ), the proof of ( 4.12 ) is similar. Since ^ is a compact, it follows that, for 
arbitrary e > 0, there exists the following e-partition of ^ 

^(e) = {y r C^:re{l,...,K}} 

such that (J r % = % , % n % ■ = 0, diam X- < £ and pick some V r £ int %, for all r^s. 
We denote 

L N (A,e) = V eW'> / exp (j3^£ <A,(o),ff,)) d^(ff.o). 

For small enough e, we have 

(l-2A^||A||e)e A '< A ^( f7 - a » <^°) < g^W^)) (l+2AT||A||e). 
Therefore, if we denote 



exp(JV(A,i?j V (<7,C7))+/3\/^£(A / (a),<7 / ))d^ v (<7,a), 



Anton Bovier and Anton Klimovsky 



23 



we get 



Hence, 



;i -2W||A||e)£ L N (T r ,A) <L N (A,e) < (1 + 2iV||A||e) £ L N (%, A). 

r=\ r=\ 



(1 -2Ar||A||e)Ljv(^,A) <Ljv(A,e) < (1 + 2A||A||e)Lyv(^,A). 
Let e 4 +0 in ( ft. 15 ) and we arrive at 

L N (A)=L N (W,A). 



(4.15) 



That is, the existence of Ljv(A) and the representation 1 3[ ) are proved. 
(2) For all N g N, we have, by the RPC averaging property (see, e.g., [^, Theorem 5.4] or Theorem ^3 
property (EJ) below), that 



1 

N" 



-E 



logL N (W,A) ='5jv(x,^)[A+iV(A,i? JV (c7,C7))] = Xb(x,.2,A,l7). 



□ 



Proof of Theorem 1.1, In essence, the proof follows almost literally the proof of Theorem 3.1 The notable 
difference is that we apply the Gaussian comparison inequality (Proposition 2.5) in order to "compute" the 
rate function in a somewhat more explicit way. 

Due to (4.5), we can without loss of generality suppose that Y is compact. For any 5 > and U g ~f, 
by ( pOC| ), there exists A(U,S) g Sym(rf), x(t/,5) g =2' (1,1) and g £'{U,d) such that 

-<17,A(I7)>+7(A(I7))<7|(17). (4.16) 

For any U g T^, there exists an open neighbourhood Y(U) C Sym(af) of t/ such that 

sup (V-U,A{U))<S. 

ve'f{u) 

Fix some e > 0. Without loss of generality, we can suppose that all the neighbourhoods satisfy additionally 
the condition diam^ / (L ? ) < £. By compactness, the covering \Jy e y Y(U) D "V has a finite subcovering, 
say Uf=i Y(UW) D "V . We denote the corresponding to this covering approximants in ( 4.16 ) by {ptS r > g 
J2'(l,l)}f =1 and {^W g £'{u( r \d)} p r=l . We have 



IlogP w (^)<Ilog(lJ/V(^(C/W)) 



Due to the concentration of measure Proposition 2.3, we can apply Lemma 3.1 and get 
1 



lim — E 



logflv (\jnu ir) ) - max E \ogP N (Y{U^)) 

\ , 7 rG{l,-,p} L 



0. 



In fact, since we know that ( 4. 1 1| ) exists, ( [4-. 18 ) implies that 



lim -E 

N]+°°N 



log/v(|J^(^ W 

r=i 



= max lim — E 

r£{l,...,p}N1+°°N 



(4.17) 



(4.18) 



(4.19) 



For t/M, * = *M, <S = Proposition |Z5| gives 



1e 

AT 



log/V(r(£/W)) 



Ie 

AT 



log JV(r(l7«))J «r)[B] 
+ ^(xW , , i/W , y (t/ (r) )) + (e) 



< -E 
~ N 



log/Vy(£/M))l -^(^^[Bj+ft, 



(4.20) 



where K > is an absolute constant. 
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By the Chebyshev inequality and Lemma 4. 1 



we have 

P N (r(U)) <P N {VeW:(V- U,A{U)) < 8} 
<e- SN f c N ^- u - A ^dP N (V) 



%"M"»^(A(l/)). 



Thus, using (4.20) and ( 4-.16| ), we get 



lim — E 



logP N (f(U^)) 



< lim 



-{u( r >,A(U (r> ))-3>[B} + i-\ogL N (A(UW)) 



ft) 



N 



-8 + Ke 



= -(U r ,A(U r )) - *[B] +I(A(U r )) -8 + K£ 
<r s (U r )-8 + Ke. 



(4.21) 



Combining ( gJj ), ( [3.1 3| ), ( [4.2 lh , we obtain 

p(^) = lim lEpogiVW] < max lim ~E [log/VW (r) )) 

< max /|(t/ (r) ) + ^e-5 

< sup/|(y)+^e-5. 

Taking 8 { +0 and e J, +0 limits, we get 

p(f) < sup ?(£/). 

The averaging property of the RPC (see, e.g., [Q, Theorem 5.4] or property (Q) of Theorem 5.3 ) gives 

p2 n 

2 



4>,v(x^)[B] = ^^x,(||G^ 1 )||^-||G«i 



k=\ 



(4.22) 



(4.23) 



To finish the proof it remains to show that, for any fixed A € Sym(c/), we have 

7(A)=X (x,J2,A,U) 



which is assured by Lemma X. 1 



□ 



4.3. Free energy lower bound. In this subsection, we return to the notations of Section X.2. 

Lemma 4.2. For any S§ C Sym(cf) such that hit J?n int 1 ^ ^ there exists Ad with hit A ^ such that 



1 

lim-E 



exp(N(A,R N (a,a)) + Y(Ai(a),<Ji))dn N (o,a) 



>log / exp(((j3 2 f/ + A)a,a))d J Li((7) > (4.24) 

J A 



Proof. In view of ( 1.10 ), iterative application of the Jensen inequality with respect to E z pt) leads to the 
following 

E[X n+1 (x,£,A,U)] <X (x,£,A,U). 
Performing the Gaussian integration, we get 

E [X n+1 (x, £,A, U)] > log ^exp (((j5 2 U + A)a, a)) dfi(a), 

where A C £ is such that ju(A) > and {R(a, ff):ff£ A^} C 33. 
Define the following Legendre transform 



□ 



r(U)= inf 

AeSym(d) 



— (U,A) — <P[B] + 1(A) +&(x, U) 



(4.25) 
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Proof of Theorem 1.2. As it is the case with the proof of Theorem |1 . l\ this proof also follows in essence 
almost literally the proof of Theorem 3.2 The notable difference is that we apply the Gaussian comparison 
in order to "compute" the rate function in a somewhat more explicit way. 

we are in the following situation: 5(~ = Sym(c/) and X is the topology 



In notations of Theorem 3.2 



induced by any norm on Sym(d). 

Let B(U,e) be the ball (in the Hilbert-Schmidt norm) of radius e > around some arbitrary U E 
us prove at first that 

1 



Let 



Iim lim-E[log/V(B(C/,e))] >/*(£/). 

ej,+0iVT°° N 



(4.26) 



Similarly to ( |4.20| ), for any (x, J3), we have 



^logP N (B(U,e)) 



N 



logP N (B(U,e)) -<P[B}+& N {x,£,U,B(U,e)) + 0(e). 



(4.27) 



The random measure P^ satisfies the assumptions of Corollary 3.1. Indeed: 



(1) Due to representation ( 4.14 ), mapping /(•) is differentiable with respect to A. Henceforth assump- 
tion ([[]) of the corollary is also fulfilled. 

(2) Let us note at first that, thanks to Proposition 13, we have = R d . Thus, the assumption ([|) 



of Corollary 3.1 is satisfied, as is condition (3.9) 



Moreover, the assumptions of Lemma 3.3 are satisfied 



(1) The concentration of measure condition is satisfied due to Proposition |2.3. 

(2) The tail decay is obvious since the family {Pn : N E \s4} has compact support. Namely, for all 
Af € N, we have suppfV = ^ ■ Thus the measure Qn.a (cf. (3.25)) generated by Pn has the same 
support. Thus, suppgw.A = . 



(3) The non-degeneracy is assured by Lemma 4.2 



Hence, due to (4.27), arguing in the same way as in Theorem 3.2, we arrive at (4.26). Note that the N \ +°° 
limit of &n(x, J2,U,B(U,e)) exists, since in (4.27) the limits of the other two ^-dependent quantities exist 
due to [ 19 1 . The subsequent e j +0 limit of the remainder term exists due to the monotonicity. 

Finally, taking the supremum over U € "V in (4.26), we get (1.22). □ 



5. GUERRA'S COMPARISON SCHEME 

In this section, we shall apply Guerra's comparison scheme (see the recent accounts by [ [l8[ |32| ||]) to 
the SK model with multidimensional spins. However, we shall use also the ideas (and the language) of 
[Jjj]. In particular, we shall use the same local c ompa rison functi onal (4.6 ) as in the AS 2 scheme, see ([5.4]). 
The section contains the proofs of the upper ( 5.16 ) and lower ( 5.18 ) bounds on the free energy without 
Assumption 1.2. The proofs use the GREM-like Gaussian processes, RPCs as in the AS 2 scheme. We 



also obtain an analytic representation of the remainder term (which is an artifact of this scheme) using the 
properties of the Bolthausen-Sznitman coalescent. 

5.1. Multidimensional Guerra's scheme. Let | = ^ (jci, . . . ,x n ) be an RPC process. Theorem 5.3 of [||] 
guarantees that there exists a rearrangement % = {^(0}/gn of the |'s atoms in a decreasing order. Recall 
( 1.16| ) and define a (random) limiting ultrametric overlap : N 2 — > [0;«] n Z as follows 



q L {i,j) = 1 +max{£ G [0;«] HZ : [&(/)]* = frUM, 



(5.1) 



where we use the convention that max0 = 0. This overlap valuation induces a sequence of random par- 
titions of N into equivalence classes. Namely, given a^eNfl [0;«], we define, for any i,j £ N, the 
Bolthausen-Sznitman equivalence relation as follows 



l 7 J 



dcf 



qh(ij) > k. 



(5.2) 
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Given n G N, assume that x and J3 satisfy ( |1.8[) and (1.7), respectively. Recall the definitions of the Gaussian 
processes X and A which satisfy (1.1) and (1.17), respectively. We consider, for t G [0;1], the following 
interpolating Hamiltonian on the configuration space Ln x srf 



H,(a,a) = VtX(a) + VT :r tA(a 1 a). (5.3) 
Given °i/ C Sym + (c/), the Hamiltonian ( |5.3[ ) in the usual way induces the following local free energy 

(p N (t,x,Q,<fr)=<S> N (x,ty)[H,}, (5.4) 



where we use the same local comparison functional ( J4.6| ) as in the AS" scheme. Using (|1 ,5|), we obtain 
then 

<p(0,x,Q,&) =<t> N (x,W)[A] and <p(l,x,Q,1f) =Q N (x,&)\X] = p N {9/). 

Now, we are going to disintegrate the Gibbs measure defined on x s# into two Gibbs measures acting 
on °i/ and stf separately. For this purpose we define the correspondent (random) local free energy on as 
follows 

exp[j3v^V// r (ff,a)j dn m (a). 



\j/(t,x,Q,a,^') = log 



(5.5) 



For a G srf ', we can define the (random) local Gibbs measure ^(t ,Q,a,^ ) G ^\ (E#) by demanding that 
the following holds 

dSf(f,jc,e,a,^; 



d\l 



(o) = lr A ,(*)(c 7 ) ex P PvNH,(o,a) - y(t,x,Q,<% ,a) 



Let us define a certain reweighting of the RPC £, with the help of (5.5). We define the random point process 
{E,}ae.e/ m the following way 

|(a) = ${a)fxp{y{t,x,Q,'%,a)). 
We also define the normalisation operation ,jV : — ► ^\(stf) as 

We introduce the /oca/ Gibbs measure ^(t,x,Q^) G x i/), for any ^ C % x &/,as follows 



^(t,x,Q,^)[r]= £ <Sr-&(a)&(t,x,Q,a,®)[y}. 



(5.6) 



Finally, we introduce, what shall call Guerra 's remainder term 
P 2 



-E 



t,x,Q,W)®&(t,x,Q,W) \\R(a l ,a^)-Q(a^,a^)\\i, 



(5.7) 



Note that (^7|) coincides with (|L20|) after substituting ( |1.18| ) with ( |L23j ), 
5.2. Local comparison. We recall for completeness the following. 

Proposition 5.1 (Ruelle [^4|, Bolthausen and Sznitman [|||). For any k G [l;n + l]flH we /iave 

E Lr(§)®^(§) |(a (1) ,a (2) ) G .fi/ 2 : </L(a\a 2 ) < it}] 

The results of Section [| can be straightforwardly generalised to the comparison scheme based on (| 
Given e, 8 > and A G Sym(c/), define 

f(A,'W,e,S) = {U'£Sym(d): \\U' - U\\ F < e, (U'-U, A) < 5}. 

We now specialise to the case = I,^(y(A,U,e, 8)). 

Lemma 5.1. We have 

^-(p N (t,x,Q,r(A,U,e,8))^(t,x,Q,L N (£f(A,U,e,8))) 
at 



(5.8) 



}2 n 



^L^(\\Q (k+l) \\ 2 F-\\Q 
L k=\ v 



(*)||2 



■0(e). 



(5.9) 
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Proof. This is an immediate consequence of Proposition 2.5 Indeed, recalling that <2(ofW, OfW) = [/, and 
setting ^ ee L(B(U,e)), we have 



2 E 



^(t,jc,e,^)®^,jc,e,^) [||^(cJ (1) ,(J (1) )-C/||^-||/?((7 (1) ,(J (2) )-e(a (1) ,a 



(2) Ml 2 



2 E 



Sf(f,jc,G,<&')®Sf(f,x ) G,<2!') ||^(cj (1) ,cJ (2) )-e(a ll, ,a w )llF 



ry(2)^|2 



-E 



^(f,^e,^)®$#(f^,e,^)[iit/i|2-ne(a( i ),a( 2 ))||2ll + ^( e ). 



(5.10) 



Using Proposition p.l|, we get 



B 2 



^( f ,x,e,^)®^( f ,x,e,^)[iit/|| 2 -iie(a( l ),a( 2 ))n 



B 2 
2 E 



?2 « 



I 



lie w iis 



?2 « 



L ||G (i+1) llS - HG W i E ® > <7l(« (1) , a (2) )} 



k=\ 



= ^Ex,(ne^)|| 2 -ii e «ii 2 



Combining ( |5TI(i| ) and ( |5.11| ), we get ( ^9| ) 
Lemma 5.2. We /jave 



(5.11) 



□ 



2 n 



/^(^(^(C/^))) ^^(^^^(BCt/^))) [A] - ^ (||<2^+D||2 _ || QW | 

L k=\ 



^(t,x,Q,Z N (B(U,e))dt + 0{e). 



(5.12) 



Remark 5.1. Note that the above lemma also holds if we substitute B(U , e) with the smaller set Y(A,U , £, 8). 
Proof. The claim follows from (5.9) by integration. □ 



Proposition 5.2. There exists C = C(E, fx) > such that, for all U G Sym + (d) as above, and all e,8 > 0, 
there exists an 8-minimal Lagrange multiplier A = A(U ', £, 5) G Sym(c/) in (1.12) smc/i that, for all t £ [0; 1], 
ana? a/Z (x, i?), we /lave 

WEjvWA^M))) < inf /(x,J2,E/,A)+C(£ + S) (5.13) 

AeSym(c/) 



and 



lim p N CL N (B(U,e))) > inf f(x,£,U,A) + lim / ,x,g,E w (B(f/,e)))df 

WT+«> AeSym(rf) ATT+~io 

-C(e + S). (5.14) 

Remark 5.2. The following upper bound also holds true. There exists C = C(E, /x) > 0, such that, for any 
A e Sym(£/), 



p N (Z N (B{U,e))) <f(x,J2,U,A) + C\\A\\ F e. (5.15) 
Pro of. The result follows from Lemma x2 by the same arguments as in the proofs of Theorems |l.l| and 



1.2 



□ 
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5.3. Free energy upper and lower bounds. Similarly to the quenched LDP bounds for the AS scheme in 



the SK model with multidimensional spins (see Section 
scheme in the same model without Assumption 1.2 on Q. 



we get the quenched LDP bounds for Guerra's 



Recall the definition of the local Parisi functional / ( |1.12[ ). 
Theorem 5.1. For any closed set V C Sym(<f ), we have 

p(Y)< sup inf f(x,£,A,U), 



where the infimum runs over all x satisfying (1.8), all =2 satisfying (1.7) and all A G Sym(af). 
Proof. The proof is identical to the one of Theorem 



1.1 



Define the local limiting Guerra remainder term 3^{x 1 J2,U) as follows 

&(x,£,U) = -1im lim [ &(t,I. N (B(U,e)))dt < 0. 

Ei+OATf+ooJO 



(5.16) 



□ 



(5.17) 



The existence of the limits in (5.17) is proved similar to the case of the AS" scheme, see the proof of 
Theorem 1.2. 



Theorem 5.2. For any open set V C Sym(c/), we have 



p(y)> sup inf \f(x,£,A,U)+@{x,£,U)], 



(5.18) 



where the infimum runs over allx satisfying (1.8); all J2 satisfying (1.7) and all A G Sym(ii). 



Proof. The proof is identical to the one of Theorem |1 .2\ The only new ingredient is Lemma |5 . 1| needed to 
recover Guerra's remainder term (5.7). □ 

5.4. The filtered <f -dimensional GREM. Given U G Sym+(c/) non-negative definite, denote by £!(U,d) 
the set of all cadlag (right continuous with left limits) Sym + (c/)-valued non-decreasing paths which end in 
matrix U, i.e., 

£{U,d) = {p: [0;l]->Sym + (d) | p(0) = 0;p(l) = U;p(t) X p(s), for t < s;p is cadlag}. (5.19) 
Define the natural inverse p -1 : Imp — > [0; 1] as 

p" 1 02) = inf{/ G [0;1] | p(t)hQ}, 

where Imp =p([0;l]). Letx = p _1 op g ^(1,1). 

Let also J2'(U,d) C J2(U,d) be the space of all piece-wise constant paths in J2(U,d) with finite (but 
arbitrary) number of jumps with an additional requirement that they have a jump at x = 1 . Given some 
p G J3'{U ,d), we enumerate its jumps and define the finite collection of matrices {2 }/^o = ^ m P C R d . 
This implies that there exist {x^lt^ C R such that 

= x < x\ < ... <x n <x n+ i = I, 

o = e< 0) * * x • • • h = u, 

where p{x^) = QW. Let us associate to p G J2'(U,d) a new path p G J2(U,d) which is obtained by the 
linear interpolation of the path p . Namely, let 



p(t) = Q^ + (Q 



(*+i)_ 



,t G [x k ;x k+ i). 



Xk+i -Xk 

Let g : R d — > IR be a function satisfying Assumption |5.l[ Let us introduce the filtered d-dimensional GREM 
process W . Let 

W = {{W k (t, [a]k)}teR + ■ a G s/,k€ [0;n] n N} 
be the collection of independent (for different a and k) IR^-valued correlated Brownian motions satisfying 



W k (?,[a] k )~(Q 



(M) _Q(k)y/2 w 



' Xk 



Xk+1 -x k 
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where {W(t)} t eR + is the standard (uncorrelated) R^-valued Brownian motion. Now, for k £ [0;n] n N, we 
define the Revalued process {Y(t,a) \ a £ &/ ,t £ [0; 1]} by 

n 

Y(t,a) = Ax k+l ,[a] k ). 



k=0 



Lemma 5.3. For a^\a^ £ we have 



Cov 



Y(t h a^),Y(t 2 ,a 



(2)< 



P [hAt2Ax qL{aW>a(2)) 



Proof. The proof is straightforward. 



□ 



Assumption 5.1. Suppose that the function g : R d — ► R satisfies g £ (R d ) and, for any c > 0, we have 
J Rd exp (g(y) - c||y|||) dy < °° and also 

sup (||Vs(y)|| 2 + ||V 2 s(y)|| 2 ) < +oo, (5.20) 

yeR d 

where V 2 g(y) denotes the matrix of second derivatives of the function g aty £ R . 

Assume g satisfies the above assumption. Let f = f p : [0; 1] x R" — > R be the function satisfying the 
following (backward) recursive definition 



f r ty) = \s(y), t = l, 

' \MogE[exp{x k f(x k+l ,y + Y(x k+u a)-Y(t,a))}], t £ [x k ;x k+1 ), 

where k £ [0;«] n N, a £ is arbitrary and fixed. 



(5.21) 



Remark 5.3. It is easy to recognise that the definition of f is a continuous "algorithmisation" of (1.11). 
Namely, X k (x,J2,A,U) = f(x k ,0), where 



f(l,y)=g(y)=log / exp(V2/3(y,(7) + (Ac7,CT}Wc7) 



(5.22) 



5.5. A computation of the remainder term. Recall the equivalence relation (5.2). In words, the equiv 



alence i ~ j means that the atoms of the RPC | with ranks i and j have the same ancestors up to the 



k-th generation. Varying the k in (5.2), we get a family of equivalences on N which possesses important 
Markovian properties, see [0], 

Lemma 5.4. For all k £ [0;n — 1] n N, we have 



.k+r 



- Xk+l ~ x ki 



(5.23) 



and also 



L^(§)(0 2 



1 —x„ 



(5.24) 



Proof. (1) To prove (5.23) we notice that 



= E 



L^)(o^)(/)-i^)(o^(oa) 



: — X k , 



where the last equality is due to Proposition 5.1 
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(2) Similarly, (5.24) follows from the following observation 



= E 



hj 
1 -x n , 



I^(§)(i)^)a)-i^(«)(o^)w 



where the last equality is due to Proposition 5.1 



□ 



Note that, using the above notations, we readily have 

1/2 N 



where {K« = {F0(i >o )} aa ,}£ j are i.i.d. copies of {F(l , a)} a e.e/- Consider the following weights 

|M(a)^^(a)exp(/(f,y(f,a))). 

— ► srf such that, for all i € INI, the following 



As in [jg], the above weights induce the permutation ft" 
holds 



|W(irW(i))>|W(irW(i + l)). 



(5.25) 



In what follows, we shall use the short-hand notations |"(?) = |"(^"(i))> (M) = ^(s, an d 

6W = {fiW(i > j) = G(* w (0»* (,) O'))}y€K. 

Theorem 5.3. Given a discrete order parameter i£ £2! ( 1 , 1 ), we /ia ve 

(1) Independence #1. The normalised RFC point process JV{£,) is independent from the correspond- 
ing randomised limiting GREM overlaps q. 

(2) Independence #2. The reordered filtered limiting GREM Y is independent from the corresponding 
reordered weights 

(3) The reordering change of measure. Given I <£ N, let Vi(-\Q) be the joint distribution of{Y(l,i)}i € j, 
and Vi(-\Q) be the joint distribution of {Y^(l,i)}i e i both conditional on Q. Then 

dv/Olfi) 



A„(\n\~n 11 e *P( x k{f{xk+i,Y(x k+l ,i))-f k (x k Y(x k ,i))}), 



(5.26) 



where the innermost product in the previous formula is taken over all equivalence classes on the 
index set I induced by the equivalence ~, 

(4) The averaging property. For all s,t € [0; 1], we have 

({« < ' > ( a '}^-e ( ")~({« ( ' , <«)}^'fi" 1 )' <" 7) 

Proof. The proof is the same as in the case of the one-dimensional SK model, see gj. □ 
Keeping in mind ( 5.26[ ), we define, for e [0;n — 1] n N, the following random variables 

T k (a) = exp (x k [f{x k+ i,Y(x k+ i , a)) - f(x k ,Y(x k , a))}). 

Given e [1;«] n N, assume that a' 1 ', a' 2 ' e .e/ satisfy ^ L (aC', a' 2 )) = fc. We introduce, for notational 
convenience, the (random) measure jx k (t , ^) - an element of (Ejy) - by demanding the following 



K(t,V) [g] = E [T\(a 1 )- ■ ■ T k (a l )T k+l (a l )T k+l (a 2 ) ■ ■ ■ T u (a l )T n {a 2 ) 

y{t,a w ,W)®y(t,a {2 \W)[g} 



(5.28) 



where g : ty 2 — > IR is an arbitrary measurable function such that ( 5.28 ) is finite. Using this notation, we can 
state the following lemma. 
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Lemma 5.5. For any i,j £ K satisfying i ~ /', i oo f we have 

k k+\ 
-1 ^2 



E[&(tJ,W)®&(tJ,W)[\\R(<j\o 2 )-Q(iJ)\\ 2 }] =pL k {t,W) \\R(a\a 2 )-Q^\\ 2 



(5.29) 



Proof. This is a direct consequence of (5.26) and the fact that under the assumptions of the theorem 

Q(iJ) = Q [k] - 

□ 

Remark 5.4. It is obvious from the previous theorem that pi k is a probability measure. 

The main result of this subsection is an "analytic projection " of t he probabilistic RPC representation 
which integrates out the dependence on the RPC. Comparing to ( 1.2C ), it has a more analytic flavor which 
will be exploited in the remainder estimates (Section 0). This is also a drawback in some sense, since the 
initial beauty of the RPCs is lost. 



Theorem 5.4. In the case of Guerra's interpolation (1.23), we have 



n-l 



)(%2 



a(t,x,Q,T. N (B(U,e))) =- £ (x k+l -x k )nt(t,I. N (B{U,E))) \\R(a\a 2 ) -Q 

L k=0 



+ 0(e) + 0{]-x n ), 



(5.30) 



as £ s-0 and x„ — > 1. 



Proof. Recalling (5.7) and (5.6), we write 

Sl(t,x,Q,T.(U,e)) = ^E[L^(|)(0^(l)W 

hj 



x&(t,x,Q,i,<&)®&(t,x,Q,j,&) \\R(a\a 2 )-Q(i,j)\\ z F 



Using Theorem 5.3, we arrive to 



L hi 



xE[^(t,x,Q,i,^)0^(t,x,Q,j,^)[\\R(a\a 2 )-Q(i,j)\\ 2 ]]. 

(We can interchange th e su mmation and expectation since all summands are non-negative.) The averaging 
property (see Theorem 5.3) then gives 

<X{t,mi,e)) = ^E[JK)(i)JK)(j)] E [&(t,i,W)®&(t,j,®) [\\R(a\a 2 ) -Q(i,j)\\$]]. 

L U 

(5.31) 

For each k e [l;n — 1] n N, we fix any indexes i'oj'o ?-/o e ^ sucn tnat 'o ~ Jo an< ^ 4 00 Jo* • ^ ear " 



ranging the terms in ( |5.31| ), we get 



\2 n 



(k) 



k=\ 



^E[^(f,/o,^)®^(f,/o,^)[||^(cT 1 ,cT 2 )-C/i]]£E[^(^)(/) 2 ]. (5.32) 



Finally, applying Lemmata 5^4 and 5L5 to ( 5.32 ), we arrive at ( |5.3C ). 



□ 
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6. The Parisi functional in terms of differential equations 



In this section, we study the properties of the multidimensional Parisi functional. We derive the multi- 
dimensional version of the Parisi PDE. This allows to represent the Parisi functional as a solution of a PDE 
evaluated at the origin. We also obtain a variational representation of the Parisi functional in terms of a 
HJB equation for a linear problem of diffusion control. As a by-product, we arrive at the strict convexity 
of the Parisi functional in 1-D which settles a problem of uniqueness of the optimal Parisi order parameter 
posed by [MM. 



Lemma 6.1. Consider the function B : R d x R + — > R defined as 

B(y,f) = ^logE[exp{x/(y + z(f))}], 



where f : R d — > R satisfies Assumption 5.1 and {z(t )}te[0;l] i* a Gaussian R -valued process with Cov [z(t)] = 
Q(t) € Sym(of) such that Q{t) u . v is differentiable, for all u,v. Then 

d t B(y,t) = i £ QuA*) $&B(y,t) +xd yu B(y,t)d yv B(y,t)) , (f,y) e (0; 1) x R d . (6.1) 

H,V=1 

In particular, the function B is differentiable with respect to the t -variable on (0; 1) and C 2 (R d ) with respect 
to the y-variable. 



Proof. Denote Z = E 



jcf(y+z{t)) 



By [0, Lemma A.l], we have 



d,B(y,t) 



1 / 1 



2x \ Z 

A straightforward calculation then gives 
1 / 1 



-E 



H,V=1 



2 P xf(z)\ 



\z=y+z(t 



d t B{y,t) 
We also have 



-E 



2x \ Z 



I Qu,(t) (Ad z J(z)d Zv f(z) +xdl Zv f(z)) e^W | 



z=y+z{t) 



u.v=l 



d yu B(y,t) = -E \x^d z J{z)\ z=y+z[t) 



(6.2) 



(6.3) 



and 



4(z" E K (Z) {x 2 d z J{z)d z J{z) + dl z J{z)) \ z=y+z{t) 



^d Zu f(z)\ z=y+z{t) ] E [xe xf Ud Zv f(z)\ z=y+z{t) ] ) 



Combining (BJ), (bj) and (6.4), we get (5.1 ). 



(6.4) 

□ 



Proposition 6.1. Denote D = UaLoC^'-^+i)- The function f = fp defined in ( 5.21 ) satisfies the final-value 
problem for the controlled semi-linear parabolic Parisi-type PDE 

'd,fM + ±L£ V=1 iPuAt) (dl yv f(y,t) + x{t)d y J{y,t)d y J(y,t)) = 0, (t,y) eDx R d , 

f(hy)=g(y), yeK d , (6.5) 

J(y,x k - 0) = f(y,x k + Q), k e [1;»] n N, ye R d . 

Note that ±p(t) = ^^- Jort 6 



Proof. A successive application of Lemma 6.1 to (5.21) on the intervals D starting from (x n ;l) gives 

Tl. □ 



Remark 6.1. Note that a straightforward inspection of (5.21 ), using ( 5.2), ( 5.3) and ( |6.4[ ), s/zows f/;af f/ze 
function f defined in (5.21) is C 1 (D) fl C([0; 1]) vv/f/z respect to the t-variable and C 2 (R d ) with respect to 
the y-variable. 
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Lemma 6.2. Given p G £2' (U ,d), the function ( |5.21[ ) satisfies the following: 

/ p (0,0) = Eflog £ |(o)exp{ g (y(l,o))}] . (6.6) 

Proof. This is an immediate consequence of the RPC averaging property ( p".27 ). □ 
Lemma 6.3. 

(1) Given k G [1 ; «] n N anc/ a non-negative definite matrix Q G Sym(af), we have 

d Q W-~ Q fp (0, 0) = - 5 (** - **-i)E [{Q,M)} , (6.7) 
where M G R rfx<: ' / 5 defined as 

M u , v =Ti(aW) ■ ■ ■ T k (aW)T k+1 (aW)T k+1 (aW)- ■ ■ T n {a^)T n (a^) 

d z u g{z)\ z=Y (\,aW) d Zv8( z )\z=Y(\,aW) 

with qi(a^', a^) = k. Moreover, ( |6.7[ ) does not depend on the choice of OC^ , OC^ G si but only 
on k. 

(2) Given a non-negative definite matrix Q G Sym(ii), we have 



d u ~ Q fp(0,0)= 1 -E[(Q,M')] 
where M' G Sym(af) is satisfies 

M; v = r, (a) • • • T„(a) (dl Zv g(z) + d Zu g(z)d Zv g(z, 



(6.8) 



as x„ — > 1. Afofe f/iaf (6.8) obviously does not depend on the choice of a G 



Proof. Applying [Q, Lemma A.l] to (5J5), we obtain 
^Eflog £ expfe(r(l,a))} 



ae,c/ 



£ ^(f)®^(f) "ft (G(a (1) ,a( 2 ')„,„| e( , )=e(i . )+se 



«,v=l 



{i a (i )=a( 2) (a^\a^)[dl ZY g(z) + d Zll g(z)d Zr g(z) 
- d zu8(z)\ z=Y (i, a W) d zAz)\ z =Y(i,a( 2 '))} 



=y(i,aW) 



j2(*)=Q« + jQ. 



Note that 



d s (e(a (1 \a (2) )«,v| G « 



'g„ lVj ?L (« (1) ,a( 2 l)=t 
0, <7 L (a (1) ,a (2) )^- 



(1) Define M(a( 1) , a (2 >) G R £/X£ ' as 

M(aW,a( 2 ))„ iV = 5 Za g(z)| z=y(1)Ce(1)) 5 Zv g(z)| z=y(liCe(2)) . 
Hence, we arrive at 

V)-q/p(°'°) = -7 E [ I l (?L(a(1 ) C( (2) )= ^(a (1) )^(a (2) )(a (1) ,a (2) )(G,M(a( 1 ),a( 2 ))) 



The proof is c oncluded si mila rly to the proof of Theorem 5.4 by using the properties of the RPC 
(Theorem 5.3 and Lemma 5.4). 



(2) The proof is the same as in (1). 

The following is a multidimensional version of [ j30| , Lemma 4.3]. 
Lemma 6.4. For any a G sf, we have 



□ 
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(1) 



dx k f P (0,0)\ 



Xk-l 



-E 



T l (a)---T k _ 2 (a)T k _ l (a)\ Xk=Xk _ 1 
(E[f(x k+1 ,Y(x k+l ,a))T k (a) 



x k= x k-i 

(2) LetM € Sym(d) withM UiV = d Zu f(x k ,Y(x k ,a))d Zv f(x k ,Y(x k ,a)), then 
.2 1 



f(x k ,Y(x k ,a)) 



Proof. This proof is the same as in [BOA. 



-E 



ri(a)---r t _ 2 (a)(e,M) 



□ 



We now generalise the PDE (6.5). Given a piece-wise continuous x G =2(1,1) andQe j3(U,d), consider 
the following terminal value problem 



d t f+±((Qy 2 f)+x(QVf,Vf)) = 0, (y,f)eR rf x(0,l), 



(6.9) 



We say that /gC([0;l] x R d ~ > R) is a piece-wise viscosity solution of (^9), if there exists the partition 
of the unit segment =:xo <x\ < ... < x n+ \ = 1 such that, for each k G [0,«] (IN,/: {x k ;x k+ \) x R d — > IR 
is a viscosity solution (see, e.g., [^]) of 

((G,V 2 /)+^(GV/,V/)) =0, (y,t) G R d x 

/(y,**+i +o) = /Cy,^+i - o), 
l/(y,i)=sG0- 

Proposition 6.2. For any p^.p^ G £!'(U,d), we have 

l/pU)(o,o)-/ p(2) (o,o)|<|| o 1 ||p( 1 )(0-p (2) (Olkd f , 

w/zereC = C(£) = E[||M|| f ]. 

Proof. This is an adaptation of the proof of [ pi] , Theorem 3.1] to the multidimensional case. Assume with- 
out loss of generality that the paths pW and have same jump times {xk}tXo- Denote the corresponding 
overlap matrices as {Q {lk) }" k to and {Q {2 ' k) }" k to- Grven s G [0; 1], define the new path p(s) G j2'{U,d) by 
assuming that it has the same jump times {jcaJ^o as the paths pO,p' 2 ' and defining its overlap matrices 
as Q( k \s) = sQ^ 1 '® + (1 — s)Q( 2 > k \ On the one hand, we readily have 

lp( 1 '(f)-p( 2 )( f )liFdf = i:(x i -^_ 1 )iie (1 ' <r) -e (2 ' ,:) iiF. 



On the other hand, using Lemma 5.3, we have 



C A 



l*/ pW (0,0)| < - I>-*n)||e (u) -G^IIf- 

L k=\ 

Finally, we have 

l/p(i)(0,0)-/ p (2)(0 > 0)|<jT |4/ pW (0,0)|dj. 
Combining the last three formulae, we get the theorem. 



□ 



Remark 6.2. Note that using the same argument and notations as in the previous theorem we get that, for 
any (y,t) G R d x [0;1], 



C(E) 



\f p m(y,t)-f p (2)M\ < 



\\pW(s)-pW(s)\\ F ds, 
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Remark 6.3. Note that we can associate to each p £ J2(U,d) a Sym + (d)-valued countably additive vector 
measure Vp £ ^#([0; l],Sym + (c/)) by the following standard procedure. Given [a;b) C [0; 1], define 

v p {[a;b)) = p(b)-p(a) 

and then extend the measure, e.g., to all Borell subsets of [0; 1]. 

Theorem 6.1. Given U £ Sym + (c/), we have 

(1) The set j3(U,d) is compact under the topology induced by the following norm 



|p(f)||fdf, p£<£(U,d). 



(6.10) 



(2) The functional^ 1 {U ,d) 3 p i-^/p(0,0) is Lipschitzian and can be uniquely extended by continuity 
to the whole 3(V ,d). 



Proof. ( 1 ) The topology induced by the norm ( 6 . 1 ) coincides with the topology of weak convergence 



of the above-defined vector measures. Since J2(U \d) is a bounded set, it is compact in the weak 
topology. 

(2) This is an immediate consequence of Proposition |6.2| 

□ 



In the next result, we summarise some results on the PDE (5.9) for the non-discrete parameters, cf. 



Proposition 6.1 



Theorem 6.2. 

(1) Existence. Assume that Q is in j2(U,d) and is piece-wise Assume also that x is in =2(1, 1) 
and is piece-wise continuous. Then the terminal value problem (6.9) has a unique continuous, 
piece-wise viscosity solution fg x G C([0; 1] x R''). 

(2) Monotonicity with respect to x. Assume Q £ j2(U,d). Assume also that i''',i' 2 ' £ =2(1,1) are 



such that jtW(f) < x^(t), almost everywhere for t £ [0; 1]. Let fgji) and fg x (i) be the corre- 
sponding solutions of (63). Then fg x (i) <fg x (i). 
(3) Monotonicity with respect to g. Assume gi,g2 ■ R d — > R satisfy Assumption 5.1 and also g\ < gi 
almost everywhere. Let fg i; fg 2 '■ ^ d X [0;1] — * R be the corresponding solutions of ( |6.9| ) with 
8 = g\, 8 = 82, respectively. Then f gl < f g2 . 



Proof. 



(1) Due to the assumptions, the diffus ion matrix Q(t) = p(t) in ( p.9| ) is non-negative definite. 
Applying [^[ Proposition 8] to the PDE ((^9) successively on the intervals [x^xt+i], where the 
p is continuous, gives the existence of the solutions in viscosity sense and, moreover, gives their 
continuity. Uniqueness is ensured by [ [l2| Theorem 1.1]. 
(2) By the approximation argument (cf. Theorem 5.1), it is enough to assume that x^ l \x^ £ J3'(l, 1) 



and Q £ £2'(IJ ,d). Then Proposition 5.1 gives the existence of the corresponding piece-wise clas- 
sical solutions of (|6^): v (i) ,/q . These solutions are obviously also the (unique) piece-wise 
viscosity solutions of (5^). The comparison result [^, Theorem 5] and the non-linear Feynman- 
Kac formula [||, Proposition 8 ] give then the clai m. 
(3) This can be seen either from the representation ( |6.6| ) and an approximation argument, or exactly 
as in (2) by invoking the results of 

□ 



6.1. T he Parisi functional. We consider now a specific terminal condition in the system ( p.5[ ) given in 

Given p £ j2{U,d), let f p : [0; 1] x R d — > R be the value of (t he con tinuous extension onto J3(U,d) of) 
the solution of ( |6.5| ) with the specific terminal condition given by ( 5.22 ). Following the ideas in the physical 
literature, we now define the Parisi functional ^(j3,p,A) : R + x £?'(U,d) x Sym + (d) x Sym(<f) — > R in 
as 



<?(P,p,A)=M0,0)-- 



*(0d(||pMi)-<C/,A>- 



(6.11) 



The integral in (6.1 1 ) is understood in the usual Lebesgue-Stiltjes sense. 
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Remark 6.4. Note that the path integral term in (5.1 1) equals /(0,0), where f(t,y) is the solution of (6.9) 
with the following boundary condition 



g(y)=p<y,i)=PZy»' y eRd - 



</ 

c 

u=\ 



Obviously £2'{d) is dense in £}(d). 
Theorem 6.3. We have 



p{P)< sup inf &>{P,p,A). 

UeSym+(d)P^'( u > d ) 

AGSym((/) 



Proof. The bound (6.12) is a straightforward consequence of Theorem 5.1 



(6.12) 



□ 



6.2. On strict convexity of the Parisi functional and its variational representation. In this subsection, 
we derive a variational representation for Parisi's functional. As a consequence, for d=\, we prove that the 
functional is strictly convex with respect to the x € =2(1, 1), if the terminal condition g (cf. (5^9)) is strictly 
convex and increasing. This result is related to the problem of strict convexity of the Parisi functional in 
the case of the SK model. 

Let W = {W(s)} se ^ + be the standard IR^-valued Brownian motion and let {J^t}teR+ be the correspon- 
dent filtration. Define 

^[t;T] = {u : [t;T] —>R d \ u is {J 5 t}igir + progressively measurable }. 

Given «6*[f;l],ge J2(U,d) and x € =2(1 , 1), consider the following Revalued and adapted to {& t } (fER + 
diffusion 

y(e,w,v)(,) =y _ J\ x ( s )Q( s )yf 2 u ( s )ds + J* (Q(s)) l/2 dW(s), »6[t;l]. 
Given some function g : R d — > R satisfying Assumption |5.l[ define /g,. x : R d x [0; 1] — > R as 



fQ,x{y,t)= sup E 



1 f 1 



2 J, 



\u(s)gds 



(6.13) 



Proposition 6.3. Let d = 1. Tfg is strictly convex and increasing, then the functional =2(1, 1) 3ih /q , 
/s strictly convex. 



Proof. We have 



y <.'*»''» , i . ( x (, ) e(,)) l / 2 M (,) di+ ^ 1 (eM) 1/2 wM. 



By an approximation argument, it is enough to prove the strict convexity for the continuous x\ ,x 2 £ =2(1, 1) 
(x\ 7^ xi). For any y e (0; 1), we have 

F (fi,Pi+(i-r)«»j) (1 ) = _y ! (^ + (1 - r )x 2 e(i)) 1/2 ^)dA + y (<2(*)) 1/2 W(*) 
< -r/ 1 {x 1 Q(s)) 1/2 u{s)6s-(l-Y)J^ {x 2 Q{s)) ll2 u{s)&s 
+ j\Q(s)) l/2 W(s) 

= r y(Q,x l ,u.t.y)^ + _ yJyteA.KW)^^ (6 14) 

where the strict inequality above is due to the strict concavity of t he sq uare root func tion. The strict 
convexity and monotonicity of g combined with the representation ( 6.14[ ) implies that ( 5.13 ) is strictly 
convex as a function of x, since a supremum of a family of convex functions is convex. □ 
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Proposition 6.4. Given a piece-wise continuous x £ =2(1, 1) and a Q G j2(U,d) which is piece-wise in 
C'(0; 1), the function fg x : R d x [0; 1] — > R defined by (6.13) is a unique, continuous, piece-wise viscosity 
solution of the following terminal value problem 

Uf+\({Q,V 2 f)+x{QVf,Vf)) = 0, (y,f)eR rf x(0,l), 

\f(yA) = g(y)- 



Proof. In a way similar to the proof of Theorem 6.2 , we successively use [ 112L Theorem 2. 1 ] on the intervals 
(x^Xk+i), where the data of the PDE are continuous. □ 



Theorem 6.4. Assume d = 1. Suppose also that g satisfies the assumptions of Proposition 5.3. For any u £ 
R, the generalised Parisi functional given by (6.11) with f p (0,0) corresponding to the terminal condition 
g is strictly convex on Q(u, 1 ). Consequently, there exists a unique optimising order parameter. 

Proof. In 1-D, we can choose the coordinates such that Q = Ut, on [0; 1]. Consequently, Q = U = const 
on [0; 1]. Hence, it is enou gh che ck the strict convexity with respect to x e =2(1, 1). The result follows by 
approximation in the norm ( |6. 10 ) of an arbitrary pair of different elements of J2(U,d) by a pair of elements 
of £>'(U,d) and Propositions jS-lj, |3]and |3| □ 



Remark 6.5. Due to the monotonicity assum ption on g, Theorem 5.4 does not cover the case of the SK 
model, where the terminal value g is given by (5.22). 



6.3. Simultaneous diagonalisation scenario. In the setups with highly symmetric state spaces (such 
as the spherical spin models of [ ^3| ] or the Gaussian spin models, see Section || below), less complex order 
parameter spaces as Q(U,d) suffice. 

Given some orthogonal matrix O E 0(d), we briefly discuss the case p 6 £!diag(U,0,d), where 

£diz g {U,0,d) = {pe £?(U,d) I for all t e [0; 1], the matrix Op{t)0* is diagonal}. 

The space i?diag(^ \0,d) is obviously isomorphic to the space of "paths" with the non-decreasing coordi- 
nate functions in R d , starting from the origin and ending at u, i.e., 

' c/ |p(0)=0;p(l) 



£(u,d) = {p : [0;1] 
where u = OUO* € 



R 



u;p(t) ^< p(s), forf < s\p is cadlag}, 



. The isomorphism is then given by 

J(M) 3 P >-> OpO* G £ diag {U,0,d). 

7. Remainder estimates 



(6.15) 



In this section, we partially extend Talagrand's remainder estimates to the multidimensional setting. 
Due to Proposition 5.2, to prove the validity of Parisi's formula it is enough to show that all the ji^ terms 
in (5.3C) almost vanish for the almost optimal parameters of the optimisation problem in (5.16). This can 
be done if the free energy of two coupled replicas of the system (7.3) is strictly smaller than twice the 
free energy of the uncoupled single system (5.4), see inequality (7.2). However, the systems involved in 
(7.2) are effectively at least as complex as the SK model itself. In Section 7.2, we again apply Guerra's 
scheme to obtain the upper bounds on (7.3) in terms of the free energy of the corresponding comparison 



GREM-inspired model. One might then hope that by a careful choice of the comparison model one can 



prove inequality (7.2). In Sections 7.3 



and 7.4, we formulate some conditions on the comparison system 



which would suffice to get inequality (7.2), giving, hence, the conditional proof of the Parisi formula, see 
Theorem [7. 1[ 

7.1. A sufficient condition for /4-terms to vanish. In this subsection, we are going to establish a suf- 
ficient condition for the measures jj,^ to vanish. This condition states roughly the following. Whenever 
the free energy of a certain replicated system uniformly in N strictly less then twice the free energy of the 
single system, the measure jU^ vanishes in — > +°° limit (see Lemma 7.2). 



Keeping in mind the definition of (cf. (5.28 )) and of the Hamiltonian H t ((j,(x) (cf. (5.3 )), we define, 
for OfW, a' 2 ' G £/^' k , the corresponding replicated Hamiltonian as 



(2) 



a 



W,a (2) ] 



(2) a (2) 



(7.1) 
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Remark 7.1. We note here that the distribution of the Hamiltonian H t (k, <J^\ a' 2 ') depends only on k and 
not on the choice of the indices a' 1 ' , a' 2 ' <G s^^< k . 



Remark 7.2. The superscript (2) in (7.1 ) (and in what follows) indicates that the quantity is related to the 
twice replicated objects. 

Define 

ee {(«« a< 2 >) g ^ 2 : <7 L (a (1) ,a (2) ) = it}. 
Additionally, for any "V C E(B(t/,e)) 2 and any suitable Gaussian process, 

{F(ct« cr (2) ,a (1) ,a (2) ) : a^ l \a^ e W\a (2) e 

we define the /oca/ remainder comparison functional as 

^ I E [io g J^IJ^ v {l3VNF(a^\a^,aW,a^)} 



(2). 



Define 



$\k,t,x,Q,r)=*®* 



(2) 



Lemma 7.1. Recalling the definition (5.4), for any Y C E(Z?(I/,e)) , we /iave 

(p£\k,t,x,Q,V) < (pg\k,t,x,Q,Z(B(U,e)) 2 ) = 2<p N {t,x,Q,Z(B(U,e))). 



(7.2) 
(7.3) 

(7.4) 



Proof. The first inequality in ( |7.4[ ) is obvious, since the expression under the integral in (7.2) is positive. 
The equality in (7.4) is an immediate consequence of the RPC averaging property (5.27). □ 

In what follows, we shall be looking for the sharper (in particular, strict) versions of the inequality ( f7.4| ) 
because of the following observation due to Talagrand [pOfl. 

Lemma 7.2. Fix an arbitrary "f C T.n(B(U, e)) 2 . Suppose that, for some S > 0, the following inequality 
holds 



<p%\k,t,x,Q,f) < 2f N (t,x,Q,I. N (B(U,e)))-e. 



(7.5) 



Then, for some K > 0, we have 



□ 



Proof. The proof is based on Theorem [2^ and follows the lines of Q21L Lemma 7]. 

7.2. Upper bounds on <p( 2 ': Guerra's scheme revisited. In this subsection, we shall develop a mecha- 
nism to obtain upper bounds on <p( 2 ' defined in (7.3). This will be achieved in the full analogy to Guerra's 
scheme by using a suitable Gaussian comparison system. 

Given U 6 Sym + (af), we say that V £ R dxd i s an admissible mutual overlap matrix for U, if 

~U V 

V* U 



e Sym+(2J). 



(7.6) 



Furthermore, define 

f(U) = {V G R dxd : V is an admissible mutual overlap matrix for U}. 

Hereinafter without further notice we assume that it € Sym + (2rf) has the form ([7~6[), where V is some 
admissible mutual overlap matrix for U . 

Let € =2(il,2d). Let y = {y/ G [0; l]}" = i be the "jump times" of the path p. We assume that the 
"times" are increasingly ordered, i.e., 

= fo <?i < ••• <?n <hi+i = 1- 
Consider the following collection of matrices 

Q = {Q / ^0( ?/ )cSym+(2 £! f)} / "+ 1 . 
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We obviously then have 



= (O) -< (1) -<...-< (n) -< (n+1) = it. 



(7.7) 



Such a path induces in the usual way the "doubled" GREM overlap kernel = {0(a^ >,w ') € 
Sym + (2<i) | a il \a i2) G ,<}, defined as 

0(a< 1 ),a( 2 ))^0^(« (1) ,« (2) )). 

We also need the <i X d submatrices of the above overlap such that 

'£l\n(a^\a^) 0|i 2 (a«a (2) ) 



0(a«,a (2 )) 



0|i 2 (a( 1 ),a( 2 ))* 0| 22 (a (1) ,a (2) ) 



(7.8) 



Remark 7.3. For (yMffP) G E^v, we sta// use f/ze notation cjW ii a' 2 ' G (IR 2 ^)^ to denote the vector 
obtained by the following concatenation of the vectors O^ 1 ' and (T^ 



i=i 



Let us observe that the process 

x (2) = | X (2) (T) =X ( CT (i)) +jr( CT (i)) | t = ffW n cj( 2 W (2) g E w } 

is actually an instance of the 2<f -dimensional Gaussian process defined in ( |l.l[ ). Hence, it has the following 
correlation structure, for T 1 , T 2 G E^, 



Cov 



\*P>(?W)\\l 



The path p induces also the following two new (independent of everything before) comparison process 
y(2) = jy( 2 )(oj) £ |R 2J | a g =04 j, with the following correlation structures 

Cov [i- (2) (a (1) ),F (2) (a (2) )l = 0(a (1) , a (2) ) G Sym+(d). 



As usual, let {Yf }f =i be the independent copies ofy< 2 ). For the purposes of new Guerra's scheme we 
define a GREM-like process (cf. ( |1.17P ) 

A( 2 ) = {A^(r,a) : t= ct^ 1 ' ii ffWjffW.fff 2 ) G E w ;a G 



A( 2 )(r,a) 



2 \ 1/2 w ™ 

I(i-/ 2) («),t ( -; 

1=1 



N 



We fix some t G [0; 1]. We would now like to apply Guerra's scheme to the comparison functional (7.2) and 
the following two processes 

These two processes are, respectively, the counterparts of the processes X(a) and A(<j,a) in Guerra's 
scheme. 

Consider a path Q G J2'(U,d) with the following jumps 

o= : e (0) ^e (1) ^...^G (n) ^e (n+1) - 

Let A = |a((7, a) : a G Ev; a G .e/n j be a Gaussian process (independent of all random objects around) 
with the following covariance structure 

E[A(cj (1) ,a (1) )A((j (2) ,a (2) )j =2(/?(a (1) ,ff (2) ),2(a (1) ,a (2) )). 

For notational convenience, we introduce also the following process 

A( 2 \a^n^ 2 Ka^\a^)^A(a^,a^)+A(^ 2 \a^). (7.9) 
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Recalling the replicated Hamiltonian (7.1) and following Guerra's scheme, we introduce, for s 6 [0; 1], the 
following interpolating Hamiltonian 



+x/r^A(v (i) Mtf (2) y i) y 2 >). 



(7.10) 



Given e, 8 > and £ € Sym(2rf), define (cf. (Q) 

r (2) (£,i!,£,<5) = {It' g Sym+(2d) : ||il'-iX|| F < e, (it'-il,£) < 8}. 
We consider the following set of the local configurations 

Z$(£,il,e,8) = {(a^a®) eL N x L N :R^(a^ n a^ 2 \a^ n a* 2 ') G r^(£,il,e,5)} . (7.11) 
Note that 4/ 2) (£,il,£,<5) C I. N (B(U,e)) 2 . We consider also the RPC £ = C(?) 

generated by the vector j: 

and, for any suitable Gaussian process 

F = {F(oW,o®,aW, a e>) | oW,o® G E^o^.oP) G <}, 
define the corresponding local comparison functional (cf. (|7~2|)) as follows 

dAi 8 V (1) )cUi V (2) )dC(« (1) )dC(a (2) ; 



Define the corresponding local free energy-like quantity as (cf. ( 5.4)) 



X{s,t,k, h Q,Q,I${Z,^e,5))=& (2 , 



>(2),*,J 



(7.12) 



To lighten the notation, we indicate hereinafter only the dependence of % on Denote 



B 



- /=] 



>('+!)! 



Lemma 7.3. There exists C = C(E) > such that, for any iX as above, we have 



Z(s,t,k,t&Q,Y§\z,il,e,8))<-&' u + Ce 



ds 



(7.13) 



Consequently, 



cp® (k, t , x, Q, 4 2) (£, ii, e, 5) ) <ogf 



E^(£,il,e,5) 



^A( 2 )(tT«nC7( 2 ),a«) 



+ VT37A(2) (cT (i) ll(T (2), a (i) ;a (2) ) 



Ce. 



(7.14) 



Proof. The idea is the same as in the proof of Theor em p.l| and is based on Proposition [2.5[ Since we are 
consid eri ng the localised free energy-like quantities ( 7.12 ), the variance terms induced by the interpolation 
( 7.10| ) in ( 2.14 ) cancel out (up to the correction 0(e)) and we are left with the non-positive contribution of 
the covariance terms. □ 

Given £ G Sym(2rf), we consider the following stencil of the Legendre transform 

& 2)M ' £ [F] = - (£,H) - B^ + lEpog J £ 2 J exp{RVNF(oM, a®, a®, a< 2 >) 

dM® w (a (1) )dM 0Ar (cJ (2) )dC(a (1) )dC(a (2) )" . 

Definition 7.1. Let F : Sym(2rf) -> R. Given 8 > 0, we call £( 0) G Sym(2rf) 5-minimal for F, if 

F(A {0) )< inf F(A) + S. 

AGSym(2rf) 



(7.15) 
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Lemma 7.4. There exists C = C(£) > such that, for all il and £ J2' (il, 2d) as above, all £,8 > 0, there 
exists a 8 -minimal Lagrange multiplier £ = £(il, £, 5) £ Sym(2d) for (7.15) such that, for all k £ [l;n] RiN, 
all t £ [0; 1], and all (x, J2), we have 



$\k,t } x,Q,Y!fi>(£,iL,e,8))< inf <J>( 2 )<^ £ V£A (2) (<* (1) n a (2) , a™) 

£GSym(2rf) 

+ VWA( 2 )(a (1) M<7 (2) ,a (1) ,a (2) ) 
+ C{e + 8). 

Proof. The argument is the same as in the proof of Theorem |l.l[ 



(7.16) 

□ 



Consider the family of matrices = joW G Sym + (2<f) | / G [0;n+ 1] R IN j, defined as 
for/G [0;£]RN,andas 

O(0 = 



g(0 g(0 



(7.17) 



(7.18) 



for / G [i + l;n+l]nM. Additionally we define, for / 6 [0; n + 1 ] , the matrices 

0«(f)=f0+(l- OH- 
Let ZW g ^2dx2d^ f OJ . i g [0;y], be independent Gaussian vectors with 



Cov 



?(0 



2)3 2 (0( /+ "(r)-0("(O 



Given y <E R , £ G Sym(2<i), consider the random variable 



45 1 (y, 3 : ) OW,£)=log / / exp( (j>Wn a^) + (£(ffW|i a^),a w n cr^) )d M (c7W)dM(tT t2j ) 



Define recursively, for I £ [n;0] R N, the following quantities 

exp( ?/ x/3(y + z(0,A ;! j;,o(0( ? ) ) £) 



X / (2) (y^, ? ,0( f ),£)^-logE 2( ' ) 



(7.19) 



(7.20) 



Lemma 7.5. We have 

$ (2),*,jx [^ A (2)( ty (i) MCy (2) ja (i)) + Vri7A( 2 )(t7W|,c7( 2 ) ! a« J a( 2 ); 



= -(£,ll)+xf (O )P) 0W(/),£). 
Proof. This is an immediate consequence of the RPC averaging property ( [5.27 ). 



(7.21) 

□ 



Proposition 7.1. Under the conditions of Lemma 7.41 we have 



<p® (M,*,e,4 2) (£,it,e,<5))< inf f-<£,il)+X { 5 2) (O, ? ,0(O,£))-B^+C(e + 5). 

£GSym(2^) v / 



,(2), 



,(2) 



Remark 7.4. Similarly to (5.15), there exists C = C(E,jtt) > 0, si/c/i that, for any £ G Sym(2<i), 

^ 2) (*,f,x,e,4 2) Wil,e)) < -(£,il>-B ? ' Q +^ 2) (O,?,0(/),£)) +C||£|| f e. 
Proof. Immediately follows from Lemmata and ^5 . 



□ 
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7.3. Adjustment of the upper bounds on <p' 2) . Proposition 2A implies that there exists r e [l;n] n IM 
such that 



lie (r - 1) i<IIF|||<ll!2 ( '' ) i. 



(7.22) 



Assume r = k. (Other cases are similar or easier as shown for 1-D in [pOj].) We make the following tuning 
of the upper bounds of the previous subsection. Set n = n + 1. Let w 6 [x r -\/2;x r ]. Define 



■, le [0;it-i]nN, 
h = Ti{w)={w, l = k, 



x h I G [i+l;n+l]nlM. 



Let 



Q 



(') = 



Q 



in 



Q { '- 1 \ I G [fc;n + 2] n IN. 



Moreover, suppose £l = {flW}?^ satisfy 



|£2<0||2 



4|ie w B, 

11= 



4||yp 



fe [0;Jt-i]nN, 

(IIQ^IIf + IMIf)- / e [jt+i ;n + 2]nN. 
Such exists due to ( 7.22 ). Moreover, if d > 2, then it is obviously non-unique. 
Lemma 7.6. In the above setup, we have 

B^^t^-x^) (||fl«|||- ||V|||) + f> (llfi^HJ- ||fi«|||)} 



Proof. The claim is a straightforward consequence of (7.23) and (7.24). 
Define the matrix € Sym + (2c/) block-wise as 

S)("+ 1 )| 11 =j3 2 f(f/-0("+ I )| 11 ) + i 8 2 (l-f)(f/-e(")) + £| 11 
D ( " +1) |i2 = j3 2 r(y-0(" +1) |, 2 ) + £| 12 , 
S) ( " +1 )| 21 =)3 2 f(y-0(" +1) | 1 2r+£|t 2 , 
®(" +1 >| 22 = p 2 t(U -&" + % 2 ) + /3 2 (1 -t)(U- Q {n) ) + £| 22 
Furthermore, we define 

"j3 2 (t/-gW)+A 







/3 2 (t/-gM)+A 



(7.23) 



(7.24) 



□ 



Sym+(2</) 3B ( " +1 ) 

Lemma 7.7. We /zave 

X^(yX{M),Z)^l°gJJ z ^p((y,<> {1) 'i tf (2) ) + (2) ( " +1) (ct (1) n ff p) ),aW ii a 

XdjLi((J (1) )dM(ff (2) )- 

Proof. Since y„ +2 = 1, the result follows from a straightforward calculation of the Gaussian integrals in 
(fL20|) for /=«+!. □ 



Define 



A 
A 



,11 = 



U qW 



Lemma 7.8. For any y G R d , / G [0; n + 2] n H we /lave 



2Xi-i(y,x,£,U,A), I G [k;n + 2] n N, 
2X { (y,x,<£,U,A), I G [0;&— 1] n N. 
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Proof. A straightforward (decreasing) induction argument on / gives the result. Indeed: for / = n + 2, an 
inspection of (|7.19[) and (1.9) immediately yields 



^+ ) 2 (3' (1) ''3' (2 ^W I Q^)-^+i(>' (1) ,^^^A)+X, 1+1 (y( 2 ^,i2,f/,A), 
where y( l \y( 2 *> G R d . Let Z( l > be a Gaussian 2<i-dimensional vector with 



Cov 



Z«l =2J3 2 (Q(' +1 )-0W). 



Define two Gaussian <i-dimensional vectors Z^- 1 and Z^' 2 by demanding that 

z(0=z». 1 i l z«> 2 . 

Due to ( |7.17[ ) and ( |7.18| ), the vectors Z^ 1 and Z^' 2 are independent, for / G [k;n + 1]. We have Z^> 
Z(')' 2 , for / G [0;fe- 1]. Assume that / G [k;n+ 1] n IN and 

^ ( +|(3' (1) ii3' (2) ,y(vv),0,£)=X / (y( 1 ),x,^,f/,A)+X / (y( 2 ),x,^,f/,A). 
By definition (7.19), we have 



x/V^^U^^logE^ 



exp| 
exp<| 



){jc I fa(y( 1 >+2^' 1 ,Jt,.2,l7,A) 

+ X / (y( 2 )+Z«' 2 ,x,^,t/,A))}" 
= X z _i(y( 1 ),x,^,t/,A)+X / _ 1 (y( 2) ,^^,f/,A). 
By the construction and previous formula, for / = k — 1 , we have 

iiy {2) ,*,?,5,£)l^, =*f(y (1) N/ 2 U,y,5,£) 

= ^_ 1 (y( I ),x,^,f/,A)+^_ 1 (y( 2 ),x,^,C/,A) 
Finally, for / G [0;£ — 2], we recursively obtain 



= ±logE*" 
= -logE- ( "' 1 



exp 
exp 



faWiyW nyW +zV\k, h Q,Z)\ w=Xk _, 



{ X j(x l+l (yM+m\x,£,U,A) 
+ X l+1 (y^+zW\x,£,U,A))}' 
= 2X,(y^\x,£,U,A)- 



□ 



Remark 7.5. Motivated by Lemmata 7.2 and 7.8 (see also Section 7.4 ), we pose the follow ing p roblem. Is 
it true that, as in 1-D (see / jjfl| |2i|/J, there exists £} G i?'(il,2af) satisfying the assumption ( 7.24 ) such that 
the following inequality holds 



inf 

£eSym(2c/) 



<2 inf 

AeSymfd) 



(£,il)+^ 2) (0, y H,5(f),£)| H . = .^ 1 
(A,f/}+Zo(0,Jc,^,C/,A))? 



(7.25) 



Similar problems have at first been posed in /[33|/. TTie resolution of the above problem seems to require 
more detailed information on the behaviour of the Parisi functional (5.1 1 ) or, equivalently, of the solution 
of (|6.9[) as a function of Q G j3(U,d) . 
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7.4. Talagrand's a priori estimates. We start from defining a class of the almost optimal paths for the 
optimisation problem in ( 5. 12). Recall the following convenient definition from 

Definition 7.2. Given U e Sy m + (d), we shall call the triple (n,p*,A*) eNx £>' n (U,d) x R d a 0-optimiser 
of the Parisi functional (6.1 1 ), if it satisfies the following two conditions 

(7.26) 



^(/3,p*,A*)< inf &{B,p,A) + 9. 

peS'(u,d) 

AeSym(t/) 

^Q3,p*,A*) = inf 3>(fi,p,A). 

pej2'„(U,d) 
AGSym(f/) 



(7.27) 



Remark 7.6. It is obvious that for any 9 > such a 9-optimiser exists. The main convenient feature of 
this definition ( as pointed out in /|3C[/J is that n ( the number of jumps of p*) is finite and fixed. 



Recalling (5.13), we set 

^,A) W ^_ (f/)A) _^£^(|| G 

^ I, 1 ^ 



2 n 



(*+l)||2. 



k=\ 



\Q^\\^+X (x,J2,U,A)- (7.28) 



Under the following assumption (at first proposed in 1-D in pCp), we shall effectively prove that remainder 
term almost vanishes on the 9 minimisers of (pTl), see Theorem P. 1. 



Assumption 7.1. Let il e Sym + (2ii) be defined by (7.6). We fix arbitrary to £ [0; 1), £ > and 8 > 0. 

There exists K = K(ta, £, 5,il) > 0, 9(t ,e,8,H) > 0, and N = N (to,E,S € N ant/ £* e Sym(2d) w'f/z 
the following property: 

If (n,p* , A*) is a 9-optimiser, for some 9 G (0; 0(foj £, 5,H)], f/ien uniformly, for all t € [0;fo), N > Nq 
and all k £ [1 ; n] fl IN, we /jave 



<pf (fe,f,x*,e*,E( / 2) (£Ml, e ,5))<20^- i? *^(f)-^||e* W -V||| + C( e + 5). 



(7.29) 



Remark 7.7. 77ze validity of the above assumption for general a priori measures is an open problem. 
However, in the particular case of the Gaussian a priori distribution the assumption is indeed effectively 
satisfied. See Section ^and Theorem 8.1 in particular. This gives a complete proof of the Parisi formula 
for the case of Gaussian spins. 

Remark 7.8. If the bound ( 7.25| ) holds then Lemma 7.6 with w = x r -\ would imply that 

q>j?\k,t,I%\2,*,il,e,8)) <20^- A *>(f)+C(£ + S). (7.30) 
The above inequality would then be a starting point for the a priori estimates in the spirit ofTalagrand /[^/ 



which might lead to the proof of Assumption 7.1 



7.5. GronwalPs inequality and the Parisi formula. 



Theorem 7.1. Suppose Assumption |Z7J holds. 
Then we have 



lim/> w (j3)= sup inf p,A). 

A 'T+~ UeSym+(d)P^'(U,d) 
AeSym(c/) 



Proof. The proof follows the argument of | ]30t l (see also pl|]) with the adaptations to the case of multidi 



Theorem 



5.1 



implies 



mensional spins. The main ingredients are the Gronwall inequality and Lemma 7.2 
that 

lim p N (B)< sup inf &>(B,p,A). 

W T+~ C/ gS ym+(rf)Pe^W) 
AeSym(d) 

We now turn to the proof of the matching lower bound. As in the proof of Theorem 1.2, it is enough to 
show that 



lim lim (p N (l,x,Q,B(U,e)) > inf 0>(B,p,A)- 

el+ONW*- pe3'(V,d) 

AGSym(rf) 



(7.31) 



Anton Bovier and Anton Klimovsky 



45 



(1) We fix an arbitrary U E Sym + (d). Fix also some fo G [0;1). By Assumption 7.1 



we can find 

the corresponding 9(to,V, U) > with the properties listed in the assumption. We pick any e 
(0;6(to,V,U)] and let («,p*,A*) be a correspondent 0-optimiser. Note that, by definition ( |7.28p , 
we have 



and, by Definition 7.2 



^(, x *,s*,AlfU_ inf &(R,p,U,A)\<6. 

peS'(u,d) 

AeSymft/) 

(2) We denote 

A N (t) = ^^^){t)-q> N {t,x*,Q*,B{U,e)). 
Note that, due to ( |5.12| ), we obviously have 

Ajv(f) > -Ce. 

Define 

A(t) = lim A N (t). 



(7.32) 



(7.33) 



The definition (7.28) and Theorem 5.4 yield 

^A w (r) < I £ (x k+l -x k )p k \\\R N (.o<U,o®)-Q®\\l 
at l k=Q l 



Ce. 



(7.34) 



(3) Let us set D = sup CTG£ \\a\\2- We note that, for any a^',a^ 2 ' G Ejv, we have 

fl((7 (1) ,CT (2) ) G [-D 2 ;D 2 ] dxd . 
Given the constant K from ( 7.29 ), for any c > 0, we define the set 



.(2),* 



(C7,e) = |VV (2) ) G E w (B(C/,e)) 2 : p^V^) -<2 W ||£ > 2^(A w (f)+c)} . (7.35) 



It is easy to see that by compactness we can find a finite covering of Y^'' k {U, £) by the neighbour 
hoods ( |7.1 1[ ) with centres, e.g., in the corresponding set of admissible overlap matrices 



(U,e)= {/?(tr« cr (1) ) G {-D 2 ;D 2 ] dxd : (a®, a®) eZ®' k (U,e)} . 



That is, there exists M = M(e, 8) G N and the finite collections of matrices {V(0};=i C ^ (17, e) 
and {U(i) }f =l C B(U,e) nSym+(d) such that 



A/ 



(7.36) 



where 



n(0 = 



C7(0 

y*(0 17(0 



GSym+(2c/), 



and £*(0 is the corresponding 5-minimal Lagrange multiplier. 
(4) Given / G [1;M] n N, let («(;'), x*(i),Q*(i),A*(i)) be the corresponding to U(i) 0(/)-optimisers. 
Due to Lipschitzianity of the Parisi functional (Proposition 6.2) and the fact that U(i) G B(U,e) 
we can assume that «(;) = «. Using the bound (7.29) and the definition (7.35), we obtain 



<2(p N (t,x*,Q*,B(U,e))-c + C(e + 8), 



where the last inequality is again due to Lipschitzianity of the Parisi functional (Proposition p.2| ) 
which allows to approximate functional's value at (jr*(0,<2*(0,A*(0) by the value at (x*,Q*,A*) 
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paying the cost of at most Ce. Choose c > C(e + 8). Then Lemma 7.2 implies that there exists 
L = L(e, 8,c) > such that 



M*(4 2) (£*AM)) <^exp 



Therefore, the inclusion ( 7.36 ) gives 



^(Ej^C/.e)) <LMexp(^-| 
Hence, for each k £ [l;n] n N, we have 



*aK<7 (I U (2) )-G w IIf 



(fc)||2 



I + II 



(1) rr(2)- 



For all cj (2) ) € (E;v(B(t/,e)) 2 \l2 ] ' k (U,e, 8) ) , we have by definition 



|/?(a( 1 ),a( 2 ))-e«|| 2 <2/:(A w (0 + C ) 



Therefore, using Remark 5.4, we arrive to 

ll<2K(A N (t)+c). 



The bound (7.37) assures that 



N 

I<LMexp[ - — 



(5) Combining fl7.39[ ) and ( [7.40D with ( |7.38| ) and ( |7.34| ), we obtain 

TV 



— Ajv(r) < 2A - (A N (t) + c) +LMexp ( - - + C(e + 5) 



dz 



Hence, 

d 
d7 



((Ajv(0 + c) exp(-2ffr )) = rap(-2») ( — {A N (t) + c) - 2K(A N (t) + c) 

-C(e + 8) 



< sxp(-2Kt) ^— (LMexp 



Integrating the above inequality and noting that due to ( 5.12 ) | Ajv(O) | < Ce, we arrive to 



A N (t)+c<(Ce + c)exp(-2Kt) +LMexp 

+ C(e + S)(exp(-2Z*) - 1) +C(e+ 8). 



(7.37) 



(7.38) 



(7.39) 



(7.40) 



Passing consequently to the limits N | +°°, £ I +0, 8 J, +0 and finally c J, +0 in the above inequal- 
ity, we get 

limA(f)<0, for all* £ [0;t Q ]. 

£|+0 

The existence of the N | +°° limits is guaranteed by the general result of Guerra and Toninelli [[T^] . 
The limits e J. +0, 8 J, +0 exist due to monotonicity. Finally, combining the above inequality with 
(7.33), we get 



limA(r)=0, for all f e [0; fo] 

£|+0 



(7.41) 
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(6) Now, it is easy to extend the validity of (7.41) onto the whole interval [0; 1]. Indeed, due to the 
boundedness of the derivatives of p# and <p, we have, for any t E [0; 1], 

M*) < / T;Ajv(/)d* 



o df 

'o 





At 



A N (t)At 



6l 



A N (t) 



dt 



<(A N (t Q )-A N (0))^ 

<A N (t Q )+L(l-t ). (7.42) 
Passing to the N | +°° limit, applying ( 7.41[ ), and then to to — > 1 limit in (7.42), we get 

limA(0=0, for all f G r0; 11. 

el+0 

(7) In particular, the previous formula yields 

= lim A(l) = ^»^*'^*' A *)(1) - lim (p N (l,x* ,Q* ,B(U ,£)). 

£i+0 £.1+0 



Note that <Pm(1,x,Q,B(U ,£)) does not depend on the choice of x and Q. Hence, by ( 7.32 ), we 
obtain 



| Mm(p N (l,x*,Q*,B(U,e))- inf 

£i+0 pe£'(U,d) 

AeSym(d) 



,p,u,A)\<e. 



The proof of ( [7.3 1| ) is finished by noticing that the 6 can be made arbitrary small. 



□ 



8. Proof of the local Parisi formula for the SK model with multidimensional 

Gaussian spins 



In this section, we prove Theorem 1.3. The rich symmetries of the Gaussian a priori distribution allow 
rather explicit computations of the Xq terms (see ( |l.ll[ )). This allows us to prove that the analogon of 
Assumption 7.1 is satisfied, implying the Parisi formula for the local free energy (Theorem 1.3). 

Remark 8.1. 77ie case of Gaussian spins is very tractable due to the (unusually) good symmetry (i.e., the 
rotational invariance) of the Gaussian measure. Therefore, it is not surprising that in this case the calculus 
resembles the one for the spherical SK model, cf. /pil BW. 



We start from the estimates under a generic (i.e., no simultaneous diagonalisation, cf. Section 6.3) 
scenario. 

8.1. The case of positive increments. Let, fork G [0;w] fl IN, 

Ae (*) = G (*+i)_ G « 

We define, for A G Sym(d), a family of matrices |z)W G R t,x</ | as follows 
and, further, for k G [0;n] H IN, 

D (t) =C-A-2p 2 f, x i A Q (l) - (8-1) 

/=* 

We assume that the matrices A and C are such that, for all / G [\;n + 1] n N, we have 

D {1) y 0. 

We need the following two small (and surely known) technical Lemmata which exploit the symmetries of 
our Gaussian setting. We include their statements for reader's convenience. 



48 



The Aizenman-Sims-Starr and Guerra's schemes for the SK model with multidimensional spins 



Lemma 8.1. Fix some vector h £ R d and a Gaussian random vector z <G R d with Varz = C 1 6 R dxd . 
Then we have 

E ; [exp «z, /i) + (Act, a) )] = (det [c (C - A) ^ 1 ] ) V2 

xexpQ^C-A)- 1 /^)). 
Proof. This is a standard Gaussian averaging argument. 



□ 



Lemma 8.2. For a positive definite matrix AQ G Sym(c/), let z ~ .yV(0,AQ). We fix also another positive 
definite matrix D G Sym(t/) such that AQ -1 >- D l . 
Then we have 



exp ( -(D \z + h),z + h) 



= (det[D(D-A2)- 1 ]) 



in- 1 / 2 



exp[-((p-AQ)-%h)). 



Proof. This is a standard Gaussian averaging argument. See, e.g., [M for an argument in 1-D 



□ 



Now we are ready to compute the term Xo(x, i2, U, A) (see ( |l . 1 l\ f) corresponding to the a priori distri- 
bution ( 1 .26 ) in a rather explicit way. 



Lemma 8.3. We have 



Proof. (1) We start from computing the following quantity 

X n+1 =log[ exp ( £ (r« , a) + (Act, a) ) d M (a), 



detD« 



(8.2) 



where FW g R d are independent Gaussian vectors with variance 



Var 



= 2p 2 AQ('\ 



We denote 



/=0 



Lemma 8.1 



gives 



/ K „ ex P ( E ,cr) + (Act, a) j djLi (a) = (det [c (C - A) " 



1/2 



xexp( - ((C-A)~7i,/z 



(2) Next, we define, for / € [Q;n] fl N, recursively the following quantities 

X / = ilogE } ''[exp(xA +1 )]. 



Applying the Lemma 8.2 to (8.2) recursively, we obtain 



Xy EE ^([(DC)]- 1 (y(°)+/ J ) ,Y®+h)+ \ gllog 



1^1. / det£>(' +1 ^ 



2,4"U/ & \ detD(') 



(8.3) 
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Recall that we have 



X = lim -logE F ° [exp(jcXi)] 

x— >+0 x 

= E Fo [Xi] 



and note that 



([D^}-\Y m +h),Y {0) +h) = 2j5 z ([D [i) ]- i ,AQ m ) + ([D^]- l h,h) 



2/rn(i)l-i An(0)\ 



Hence, combining (8.4) and (p. 50 with (p. 3), we obtain the theorem. 



(8.4) 



(8.5) 



□ 



8.2. Simultaneous diagonalisation scenario. In what follows, we employ the simultaneous diagonalisa- 
tion scenario introduced in Section [T^. Suppose that, for / G [0; n + 1] n N, and some matrix O G 0(d), we 
have 

D {l) = *d {l) 0, 
where the vectors d {,) G R d , for / G [0;n] n N, satisfy 

That is, the vectors are (component-wise) increasingly ordered and non-negative. 
Lemma 8.4. We have 



Xb(x,J2,tf,A) = -£ 



2 

*• k=\ 



4 



It: 10 



(/) 



J 2 n 



I** (lie 



( * +1) iiJ-iifi (t) iiJ 



;2 n 



I>(l 



^lli-ll^lli 



(8.6) 
(8.7) 



Proof. This is a standard argument which relies on the standard invariance properties of the determinant 
and the matrix trace. 

□ 



Define the 1-D Parisi functional for the case (1.26) as 



2J3V 1 >+/j 2 " 1 d^ 



(J 2 E" (\i ,M) t-W"\ 

1=1 



(8.8) 



Proposition 8.1. There exists C — C(Y) > such that, for all u G ¥> d + and all e , 8 > 0, f/iere ejc/sfs an 
8-minimal Lagrange multiplier X = X(U,E,8) G R d in ( 1.12 ) such that, for all t G [0; 1] ana? all (x,p), we 
have 



p N (Z N (-r(A,U,e,8)))<-mf[ £ &(p v ,K) ) +C(e + 8) 



(8.9) 



and 



lim PN (L N (B(U, e)))>Unf I V^(p v ,K)+ Hm f 1 ^(t,x,Q,L N (B(U,e)))dt 

N]+oo 2 p.X N]+^Jo 

+ C(e + 8). 

Proof. We combine ( |8~6| ) and (Q) and the Proposition to get ( |8~9| ) and ( |8.10| ). 



(8.10) 

□ 
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8.3. The Crisanti-Sommers functional in 1-D. In this subsection, we adapt the proof of [ p9| | to obtain 
the equivalence between the (very tractable) Crisanti-Sommers functional flllj ] and the Parisi one (3.8) in 
the case of the Gaussian a priori measure (1.26). Similar ideas based on the symmetry of the a priori 
measure were exploited in the case of the spherical models by p| |23|. 

We restrict the consideration to 1-D situation for a moment. Given u > 0, consider p e 1), A 6 R, 

AeR and let {c/ (/) G R}'/!/ be the scalars playing the role of matrices dW ( c f. |y])). That is, 



d {l) =c-X-2$ 2 Y^x k L 

k=i ^ 

dW = c. 

We define, for k G [1;«] n N, the family of vectors {*(*) G R d }^ =0 by 



.(*) = 



l=k 



(8.11) 



We also define the Crisanti-Sommers functional as follows 

(1) B-l 1 

^(P) =1- CM +/ ! V 1) + i--+^ -log 



Jl) 



,(/+!) 



log 



c(u-?W) 



/=i 



Lemma 8.5. If (p,A) is cm optimiser for ( |8.8| ), nW is, 



^(p,A)= inf ^(p',A' 
(p'.A') 



then, for all k G [1;«] n N, the pair (p,A) satisfies 



(k) _ h 2 +2p 2 qW k fl I f I l 



Moreover 



A =c-2j3 2 ( M - 9 ( " ) )-(«-9 ( " ) ) _1 , 



onrf, /or all k £ [l;n] fl N, we /zave 



.(*) 



(8.12) 



(8.13) 



(8.14) 



(8.15) 



(8.16) 



(8.17) 



Remark 8.2. /« the formulation of the theorem (as well as elsewhere), it is implicit that = d^ k \p,X) 
arids® =5fW(p,A). 



Proof. (1) Rearranging the terms in (p.8|), we observe that 

1 1 



^- 1 



ll 0g j(«+i)_ll 0grf (i) 

A/2 -^1 



1=1 

We compute, for fc, / G [l;n] D N, 



o, 



(8.18) 



(8.19) 



2p 2 (x k -x k -i), k>l. 
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Using ( 8.1 9| ) and the representation (8.18), we compute the necessary condition for (q,X) satisfy 
(|L27|), for k £ [2;«]nN, 



= & (q, A ) =2/3 2 (x k - x k _ , 



i 



■qk 



We also have (for k = 1) 



0: 



^(?,A)=2j3 



(8.20) 



2/3 ^ 



(i) 



[dW] 2 

qW+h 2 



X] 



Xl i (1) 



(8.21) 



Relations Q8.20D and d8.21| ) then imply d8.14| ). 
(2) Using the fact that 
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we obtain 



-^( g ,;i) = - M + 



1 



1 



(8.22) 



Applying ( 8.14 ) with k = n in ( 8.22 ), we obtain that the necessary condition for A to satisfy ( 8.13 ) 
is as follows 

0=-£-£»(?,A) = -u + g M + _L (J___L^ 
<9A v y jc„ WW 



-u + q y 



1 



M + < ? ( " ) +fc-A-2j3 2 ( M - ? W) 



(8.23) 



which implies (8.15). 



(3) Relation (8.16) is proved as follows. Subtracting the relations (3.14), we obtain, for£ £ [1;« — 1] n 

K 

l l 



Xk (q [k+l) -q {k) 



By ( |8.23 ), we have 



(8.24) 



^(^ +1) -^)=»-^ = i 



(That is, ( |8.24 ) is valid also for A; = n.) Combining the previous two relations, we get, for k £ 
[l;«]nN, 

1 



(8.25) 



Using ( |8.25| ) and (|8_24j), we get 

2jS 2 JC,(^ +1 )- 9 «)=^ +1 )-d« 

^)) =d( k+ VdWx k U k +V-qW) =d^dW(s^-s^) 



(by 



1 1 



which is (8.16). 



□ 
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Lemma 8.6. If p is an optimiser of (8.12), that is, 



P' 



then, for all I G [1 ; n] n H ( 8. 16 ) holds. 



Proof. The strategy is the same as in the previous lemma. We rearrange the summands in (8.12) to get 



^ p =h 2 s W + ^ry + — — 

sW x\ Xn-\ 



1=2 



We have, for k, I G [1 ; n] n I 



.log(c( M -^W))+ J 8 2 £x / ([^ +1 )] 2 - 




(8.26) 



(8.27) 



(1) Relation (8.27) implies, for k G [2;n-l]nN, 



(i) 



tf,y(p) =hT{x k -i -x k ) - 77^(^-1 -x k ) + 



t ,(l)]2 

2^V* ) (*fc-i-*k)=0. 



V 1 ) 



L 

/=2 



r(*) 



x k in 



Hence, 



2 „« 



i-1 



1 / 1 



/=2 



Xi Xl-\ 



= -h l 



q 



(8.28) 



(2) To handle the case k = n, we note that 

logfl+c(«-?W) 



■log 



and, hence, the argument in the previous item shows that (8.28) is also valid for k = n. 
(3) Differentiating the representation (8.26) with respect to g' 1 ) and using (8.27), we obtain 



d ^(p) = -^ + JL.^ (1 



dqO) 



Therefore, 



2.(1) 



sW [sW] 2 XisW 
(i) 



-2j3 2 jci^ (1) =0. 



l2 , T 



(1)12 



which is (8.28), for it = 1. 



(4) Subtracting equations (8.28), we arrive to (8.16), for all k G [l;n] D N. 



□ 



Proposition 8.2. The functionals (8.12) and ( |8.8[ ) are equivalent in the following sense 

p'.V p' 



Proof. (1) Let (p,A) be the solutions of equations (8.16) and (8.15). Le mma 8.6 guarantees that 
p is the optimiser of the Crisnati-Sommers functional and Lemma SL5 assures that (p,A) is the 
optimiser of the Parisi functional. 
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(2) We have 



(i) 



, A) - tf.Y(p) = - Am + 2/3V 1 V 1 ' - * +cu _ 1 



2fl 2 L*([9<' +1 >M9°>] : 
/=i v 



(8.29) 



We can simplify the <t>[B]-like term (that is the summation) in (8.29), using ( fj6[ ) and (|T|). 
Indeed, 



n-l 



n-l 



/=i /=i 

n-l / 

(by ( |8~T6| ) and ( fIT| )) = £ 2/3 V +1) 



/=i 



9 



(0 



1 1 



(8.30) 



Regrouping the summands in ( |8.30| ), we get 

( pOl ) =2j5 2 £ a O + 2/J 2 (*< 1 V 1 > 



/=i 



/=! 



Due to (8.16), we have 



2j3 2 (V' +1 >-^ 



9 w_ 9 m 

.(») ,(i) 



,(/)_«.(/+!) „(/+!)_„(/) 



(8.31) 



Applying the previous relation, we get that the both summations in (8.31) cancel out and we end 
up with 



,(«) 



Now, turning back to ( J8.29 ), we get 

^(p,A) - Vy{p) = -Am -2/3 2 (ii 2 - [<? (n) ] 2 ) +2p 2 q { " ) s^ - g- + CM - 1 
(by ( p5b ) and ( pH ) = -m (c - 2/3 2 (m - $(»)) - (u - q^)' 1 ) - 2/3 2 (ii 2 - [q {n) \ 

- — ~~~TT + 2/3 2 t? (") f ii - q [n) ) + cu - 1 



□ 

8.4. Replica symme tric c alculations. In this subsection, we shall consider the one dimensional case of 
the a priori measure ( 1.26 ) with h = 0. We shall also restrict the computations to the case n = 1 which 
is often referred to in physical literature as the replica symmetric scenario. It is indeed the right scenario 



under the above assumptions, as shows Theorem 1.3 



Lemma 8.7. Let ji satisfy (1.26) with h = 0. Assume d = 1, n — 1 and c > 0. Given u > 0, we have 



qe[0;u] 



inf <€&{p) = inf ( 1 -cu + log ( c (u-q)) + -^—+ p z (u z - q z ) ) = f{c,u) 



?2/„2 „2i 



where f(c,u) is defined in (1.27). 
Proof. Using the definitions, we obtain 



dq dq 



u — q 



(8.32) 



\og{u-q) + — !— + j3 z (u z -q z 
u — q 



{u-q} 



-2p 2 q. 
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Hence, the critical points of q i— ► c €S fi {q, u) are 



V2 

qo = 0,qi. 2 = u±—. 



Furthermore, we also have 



1 



2q 



{u-qf {u-qf 



2j3 2 



Hence, as a simple calculation shows, the infima in (8.32) are attained on 

2/3 



which implies (8.32). 



(8.33) 

□ 



Lemma 8.8. Under the assumptions of Lemma $.7\ we have 
(1) Forc>2^/2fi, we have 



sup inf tfy(q,u) = tf<y(0,u* )= j3 1 (u*) z +log cu* - cu* +1, 

„>(HG[0;u] 



where 



(2) For c < 2V2/3, we /lave 



sup inf c <o.y{q,u) = +°°. 
„>0</g[0;h] 



Remark 8.3. Under the assumptions, the above theorem says that from the point of view of the global 
free energy, the system can only exist in the "high temperature" scenario, cf. (1.27). The threshold at 
cq = 2\/2j3 could be easily understood from the perspective of the norms of random matrices. 

Proof. (1) Suppose c > 2-^2/3. Recalling (1.27), for u G (0; we introduce the following function 

f(u) = log(cw) + J3 2 u 2 - cu + 1 . 

We have 



4-/(m) = -+2j3 z w-c. 
au u 



Hence, the critical points of the function / are 



"1,2 



c±y/c 2 -&p"< 
4p 2 



Furthermore, we have 



/( M )=2/3 2 --. 
au A u 



We notice that u* < ~& and, hence, due to ( 1.27 ) 

tfy(0,u*) =J3 2 (m*) 2 + log cu*- cu* + 1. 
(2) If c < 2a/2)3, then the function 

c 1 

«h (2V2j3-c)w + log-^--(l+log2) 
is unbounded on +°°). 



□ 
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8.5. The multidimensional Crisanti-Sommers functional. Recall the definition ( 7.28[ ). 
Proposition 8.3. Assume d = 1. Given u > 0, we have 



20 



(.Y*,j2?*,A*) 



(0 



3V2p-c) M + logj|-l-^-/(y2 M j3- 
2j3 2 (m) 2 + log(c M ) - cu + 1 - fj3 2 0) 2 , 



M > 



M < 



V2 



Proof. Combining (8.8), (8.12) with Lemma 3.7 and Proposition |8.2j, we get the claim. 



(8.34) 



□ 



8.6. Talagrand's a priori estimates. In this subsection, we prove that Assumption 7.1 is satisfied in the 
case of the Gaussian a priori distribution (1.26) with h = 0. 



Theorem 8.1. Let pL satisfy ( 1 .26 ) with h = 0, assume U G Sym + (d) is such that min v u v > ^ and suppose 
OQ.LetQ^ Q* and A = A*. 



Then, for any to G (0; 1) and any t G (0;fo], we have (cf. (7.29) with k = I) 



4 2) (l,/ J x,<2,4 2, (£,il,e,5))<20^ A )(/)-^||(2W 



42), 



1 



(8.35) 



Proof. (1) We employ the notations of Section 7.2. Let n = 1. Given il G Sym(2c/) (cf. (7.6)), 
choose arbitrary matrices G Sym(2ii) | / G [0;2] n n| satisfying ( |7.7| ). Define y = jc which, 

in particular, implies that £ = Finally, we set, for / e [0;n + 1] n N, = Q^> . 
(2) Proposition 7T implies that, for any 5-minimal £ G ^} dxld , we have 

f/3 2 



^ ) (l,t,i 1 fi ) ^ ) (£,U 1 e 1 «)) < -(£,U) ■ 



\Q {2) \\1- 



in (i) nr 



(3) We define a matrix £ G R 2c/xM as follows 



-zf (l,r,fi«(*),£) + ^(£ + S). 



C 
c 



(8.36) 



Recalling (|Q|), we define also the following matrices £)' 2 ' = £ and 

s) (1) = e:-£- (o (2) (0-£ (1) (0^ 



(8.37) 



Applying Proposition 8.1 to (B.36), we get 



4 ] (1 , t , 4 2) (£,il, e, 5)) < i [- (£, It) - r/3 2 ( ||0 (2 > \\ 2 F - 



+2)3 2 ([D( 1 )]- 1 ,5( 1 )(r))+log 
:$( 2 )'^' £ + ^(e). 



(i) nr 

/ detS< 2 )' 
\ detDf 1 ) 



■0(e) 



(4) Assume that the matrices 



Q(l) ;0 (2) ;2) (l) eR 



2dx2d 



(8.38) 



(8.39) 



are simultaneously diagonalisable in the same basis which is given by the orthogonal matrix £) G 
U 2dx2d . Let the vectors 



(8.40) 



be the corresponding spectra of the matrices (8.39). That is, we assume that 
O' 11 =D*diagq( 1 'D,0 (2) =D*diagq( 2 )£), 
= £)*o (1) D,Q (1) = D*diag0'( 1) £), 
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where we have introduced the matrix 0'W(f) e Sym+(2c/). By ( |8.33[ ), we have, - = ^> 
where / denotes the unit matrix of the suitable dimension. The definitions (7.17) and (7.18) then 
imply 

£(2)_0(l) = ^L 

2j3 

Using the definitions and the above relation, we obtain 

ffl\t) = 0*(?diagq« + (1 -t)Q'W)Q, 



(8.41) 



& 2 \t) - 0«(/) = O* (t diag(q( 2 > - q«) + (l-t)^l)D. 



Motivated by (8.17), we set 



fv = ( Uv - q,\ 

In view of (8.37), the above choice necessarily yields (cf. ( |8.15 )) 



(8.42) 



(8.43) 



£ = £- £)*diag(u„ - qi 15 )- 1 ^ - (fi (2) (0 



e:-D*(diag(u„-q! 1) )- 1 +?diag(q( 2 )-q( 1 )) + (l-f)I-/]0 



(8.44) 



Applying Lemma 3.4 to ( p.38| ) and using (3.44), ( 3.43), ( p.42| ), we get the following diagonalised 
representation of ( 3.36) 



<$ ] ( 1 , t ,x, Q, E^ 2) (£,ii, e , 8)) < l - log det £ - l - (£, it) 



V5 



2j3 z (u,.- q, (1) ) (fq, (1) + (1 -f)qv 1J ) +log(u,.-0^O 
r/3 2 ((q, (2) ) 2 -(q 1 (1, ) 2 )} + ^). 



(l) 



2/3 

Mi 



(8.45) 



Using the definitions, we get 



(£,11) = 2(C,I/)=2£c v k v 



log det £ = 2 log det C = 2 £ \ogc v 

v=l 



Motivated by (8.41) (or by (8.33)), we define 



q«:=u, 



V2 
v 2j8 ■ 



(8.46) 



(8.47) 



In this case, as a straightforward calculation shows, the expression in the curly brackets in (8.45) 
equals 



2V2j3u v + J3 V53$ } (1 - - log/5 - ^ (log 2 - r ) . 
By the definitions and the general properties of matrix trace, we have 

2d 2d d 

v=l r=l v=l 

2d d 

^Uv = 2£C7 v , v . 

1=1 v=\ 



(8.48) 



(8.49) 
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Combining ( |8.45| ) with ( |8.48| ), ( |8.49| ) and ( |8.46| ), we obtain 

d 



( P ^\l,t,x,Q,^\£,il,e,8))<'£(-CvU v +logc v + 3V2uP 

v=l 

- i(log 2 - f ) - V2>m - log/3-1 
= 2l>('fav, + tf(e), 



(£) 



(8.50) 



where in the last line we have used the relation (3.34). 
(5) To get the version of the a priori bound (8.5C ) with the quadratic correction term as stated in ( 8.35| ), 
we perturb the r.h.s of ( 8.36 ) around our choice of 2)W in ( |8.43 ), i.e., 



where in the last equality we used ( |8.47[ ). 
8.7. The local low temperature Parisi formula. 



□ 



Proof of Theorem 1.3. The result follows from Theorem 8.1 and Theorem 7.1. Note that the proof of 
Theorem 7.1 requires a minor modification to cope with the fact that the a priori distribution (1.26) is 
unbounded. This minor problem can be fixed by considering the pruned Gaussian distribution and using 
the elementary estimates to bound the tiny Gaussian tails. □ 



Appendix A. 



The general result of Guerra and Toninelli 1 19] implies that the thermodynamic limit of the local free 
energy (1.6) exists almost surely and in L . The following existence of the limiting average overlap is an 
immediate consequence of this. 

Proposition A.l. We have 

E [SfjvQS) ® Sfo(/3) [VmH n {o) - E [H N {o)H N (o')]]] — -» C(j3) > 0, 

where C : R+ — * R+. 

Proof. The free energy is a convex function of /3 (a consequence of the Holder inequality). Hence, by a 
result in Jl6] ] the following holds 

Um^E M /3)] = A Eb(/5)] . 



Proposition 2.4 implies 
d 



d/3 



E [pjv(j3)] = j3E ]efs{fi) ®^v(J3) [Var//;v(a) - E [H^(ct)« v (o / )]] 



□ 



The follo wing super-additivity result is an application of the Gaussian comparison inequalities obtained 
in Subection 2.3. Note that the result does not provide enough information for the cavity-like argument of 
§■ 

Proposition A.2. For any = B(U, e) C °i/ , we have 

NE[p N (r)]+ME\p M (y)} < (N+M)E\p N +M{y)] + {N+M)&(e), 

as e i +0. 
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Proof. Define the process Y^,m = {Y{(j) : C7 = Of n T; a £ T,n,? € £*f} as follows 



Y(a ii t) : 



N 



N + M 



1/2 



4V) + 



M 



1/2 



whereX^ andX^ 2 ^ are two independent copies of the process X. Given some Gaussian process {C(<j)} ae -£ N , 
let us introduce the functional <&#(/J)[C] as follows 

<W0)[C] = E [lo gA i®( w+M ) [l Iw(n l lM(r) exp( i 8V5vTMC) 



Now, set (p(t) = <&n+m{P) \y/iXtf + M + V 1 — tY NM \ ■ Applying Proposition 2^5, we get 



Cov 



Yn,m(0 {1) ),Yn,m(o {2) ) 



- (Cov [x N+M (a^),X N+M (a^) 
Note that we have 

<p(0) =NE\p N (y)]+ME\pn{r)], 
(p(l)<(N+M)E\p N+M (f)], 

where the last inequality is due to the fact that, for all a £ TLn{Y) and all X £ Ejv("^)> we nave 

anTeE Ar+M (' , / / ). 

Moreover, for a = a ii X with a £ and a £ Em(^) we have 



(A.1) 



(A.2) 



VarX;v +M (c7) - VdxY NM {o) 



N . . M , s 
R N (a, a) + tz r m(*, 



N + M 
M 



N + M 



N + M' 



R N (a,a)\\: 



N + M 

Also, due to convexity of the norm, we have 

X N+M {o^),X N+M {o [ 



\R M {x,x)\\l = G{e). 



Cov 



N 



N + M 

M 
N + M 



^(« (1) ,« (2) ) 



-Cov 
M 



N + M 

\Rm{x [X \x (2) )\\1<Q 



Y n ,m{<? {1] ),Ynm{o {2) ) 
,(* ( V 2) ) ' 



N 



N + M 



R N (a {1 \a 



W /y(2)M|2 



Applying / df to ( A 1 ) and using the previous two formulae, we get the claim. 



□ 
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